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CHAPTER  1 


INTRODUCTION 


1.1.  Motivation 

Many  physical  systems  are  modeled  adequately  by  a  system  of 
ordinary  and/or  difference  equations.  However,  the  need  to  consider  models 
with  partial  differential  equations  arises  in  many  areas  of  the  phvsical 
sciences.  Examples  include  mass  or  heat  transfer  [10,14,38],  eJ  cicity 
[10,33,34,37,38,40],  electromagnetic  wave  propagation  [10,37,38  luclear 
reactor  theory  [9,34],  fluid  flow  [19],  and  stochastic  processes  ^  ,. 

Due  to  physical  considerations  and  the  desire  to  obtain  simpler 
models,  the  engineer  or  the, applied  mathematician  often  redefines  the 
variables  of  the  model  at  hand,  so  that  small  parameters  appear  explicitly. 
In  systems  described  by  partial  differential  equations,  the  small  parameter 
may  represent  a  small  diffusivity  or  a  small  convection  coefficient  in 
heat  transfer,  the  thinness  of  a  vibrating  membrane  in  elasticity,  or 
"cheapness"  of  control  in  optimal  control  problems.  However,  the  introduc¬ 
tion  of  small  parameters  may  be  purely  artificial.  Such  is  the  case  in 
regularized  and  penalized  problems  [23].  In  this  situation,  the  interest 
usually  lies  in  the  properties  of  the  solution  of  the  limiting  problem  as 
€ -*■  0  and  not  the  problem  itself. 

The  dependence  of  these  models  on  c  is  singular,  i.e.,  the  formal 
limit  of  these  operators  as  may  or  may  not  exist.  In  the  case  it 

exists,  ellipticity  of  the  original  operator  is  often  lest  (order 
reduction) .  Hence  the  solution  of  a  boundary  value  problem  involving  the 
perturbed  operator  converges  to  the  corresponding  boundary  value  problem 


involving  the  formal  operator  limit  in  a  larger  Hilbert  space.  Thus,  the 
need  to  introduce  "correctors"  concentrated  in  the  vicinity  of  the  boundary 
of  the  set  over  which  the  Hilbert  space  of  functions  is  defined,  is 
inevitable  [1,5,11-14],  [16,24,25,27,29,31,41].  In  some  instances,  the 
formal  operator  limit  may  not  be  well-defined  (e.g.,  not  elliptic). 
Depending  upon  the  type  of  boundary  value  problems  at  hand,  their  solutions 

may  be  expanded  in  Laurent  series  expansion  of  c  or  may  be  approximated 
by  "regular"  expansions  with  the  addition  of  correctors  [24,27].  In  this 
class  of  problems,  if  ellipticity  of  the  formal  operator  limit  is  not 
uniformly  lost  (e.g,,  in  some  stiff  operators),  one  would  expect  some 
"separation"  in  the  spectrum  of  the  original  operator.  This  last  intuitive 
observation  is  the  driving  force  behind  the  present  investigation  of  the 
spectral  decomposition  of  stiff  operators. 


1.2.  Literature  Survey 

Asymptotic  expansions  of  linear  and  nonlinear  differential  operators 
depending  on  a  small  parameter  e  have  been  studied  by  several  scientists 
over  the  past  several  decades.  The  underlying  theory  (such  as  order  and 
validity  of  approximation,  asymptotic  error  estimates,  etc.)  is  discussed 
in  detail  in  [11].  Eigenvalue  problems  of  some  of  these  operators  are 
considered  in  [5,12,13,15,16,31,32]  and  the  references  therein.  Most  of 
these  references  assume  that  the  formal  limits  as  c  ^ 


0  of  such  operators 
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are  uniformly  elliptic.  In  this  instance,  the  eigenvalues  are  uniformly 
bounded  away  from  zero  as  e-^0,  or  may  become  dense  in  a  subset  of  the  real 
line  [13]. 

In  stiff  operators,  there  are  several  conditions  that  cause 
stiffness,  some  of  which  are  explored  in  this  thesis.  In  general,  the 
eigenvalues  (and  the  corresponding  eigenvectors)  of  stiff  operators  can  be 
decomposed  into  groups  depending  upon  their  convergence  as  e-*-0. 

Boundary  value  problems  for  some  classes  of  operators  that  depend 
upon  a  small  parameter  e  (including  several  control  problems)  are  studied 
in  [1,11,14,24,27,29],  just  to  name  a  few.  Formal  asymptotic  expansions  of 
the  solutions  of  stiff  elliptic  boundary  value  problems  are  considered  in 
[24,27,29]  without  any  reference  to  their  spectral  decomposition.  Parabolic 
and  hyperbolic  problems  involving  stiff  operators  seem  not  to  have  been 
previously  investigated. 

The  main  contributions  of  this  thesis  are: 

1)  The  spectral  analysis  of  some  stiff  operators,  including  the 
convergence  of  their  eigenvalue-eigenvector  pairs  as  e  ->•  0.  An  appropriate 
terminology  such  as  flattening,  attenuation  and  oscillation,  is  introduced 
to  describe  the  deformations  of  the  eigenvectors  as  e  0. 

2)  The  approximation  of  the  solutions  of  boundary  value  (including 
some  control)  problems  involving  stiff  operators,  using  the  weak  limits  of 
the  eigenvectors  of  the  aforementioned  operators.  The  advantage  of  this 
method  is  its  simplicity.  Moreover,  it  alleviates  stiffness  when  numerical 
algorithms  are  employed. 
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1.3.  Thesis  Overview 

In  Chapter  2,  the  eigenv.  l-'.e  problem  of  stiff  operators  that  have 
coefficients  0(1) ,0(e)  , . . . ,0(e^)  in  the  different  interfaced  subsets 

, . . .  ,fip  (whose  union  constitutes  the  open  connected  set  2  (Cr’^)), 

are  analyzed.  The  interfaces  are  the  counterpart  of  interconnections 
between  "areas"  in  large  scale  lumped  systems.  One  way  to  understand 

the  singular  behavior  of  stiff  systems  is  to  analyze  their  spectral 
decomposition.  Indeed,  for  small  values  of  e,  the  eigenvalues  of  stiff 
operators  can  be  separated,  depending  upon  their  convergence  as  c-*-0.  Their 
corresponding  eigenvectors  are  also  classified  accordingly. 

In  Chapter  3,  using  the  convergence  results  of  the  eigenvalue- 
eigenvectors  of  stiff  operators  as  e->0,  approximations  to  some  classical 
boundary  value  problems  (namely  elliptic,  parabolic,  and  hyperbolic)  are 
constructed  and  asymptotic  error  estimates  are  derived.  Most  of  the  ideas 
are  specialized  to  second  order  operators  for  simplicity.  However,  the 
approach  is  general  enough  and  hence  may  be  applied  in  many  other  similar 
problems . 

In  Chapter  4,  two  control  problems  are  considered  using  the 
approach  developed  in  Chapter  3.  For  optimization  problems  (including 
control  problems),  some  caution  is  advised  in  applying  this  approach, 
because  of  the  inherent  dependence  of  the  optimality  systems  on  the  type  of 
observation  and  control  [23,29].  There  are  several  control  and  observation 
mechanisms,  i.e.,  distributed,  boundary,  pointwise,  etc.  The  control  action 
transforms  the  characteristics  of  the  system.  Consequently,  it  may  not 
be  possible  to  use  the  eigenvectors  of  the  uncontrolled  system  to  solve 


the  controlled  one. 


Chapter  5  gives  the  numerical  results  concerning  boundary 
value  and  control  problems.  The  approximations  derived  using  the  approach 
of  Chapter  3  are  compared  with  those  of  the  direct  approach  using  a  finite 
element  method. 

The  last  chapter  contains  some  concluding  remarks  and  some 
possible  extensions  of  the  results  presented  in  this  thesis. 
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CHAPTER  2 

SPECTRAL  ANALYSIS  OF  STIFF  OPERATORS 

2.1.  Introduction 

This  chapter  considers  the  eigenvalue  problem  of  the  following  formal 
selfadjoint  operator,  written  in  matrix  form  as: 


where  A^,  i  =  0,1  are  unbounded  operators. 

Many  physical  nrobleras  can  be  described  by  models  containing  the 
operators  A^.  Examples  of  such  problems  in  distributed  param.eter  systems 
are  numerous.  Without  being  exhaustive,  examples  include  the  following: 

1)  Nuclear  reactor  operations  [9] 

2)  Keat  or  mass  transfer  in  interfaced  media  having  different 
dif fusivities  [10,38] 

3)  Electromagnetic  wave  propagation  in  waveguides  made  of  materials 
having  different  permittivities  [10] 

4)  Small  vibrations  of  elastic  interfaced  through  membranes  with  different 
material  densities  [40] 

5)  Continuous  stochastic  processes  when  the  noise  intensity  level  is 
different  from  one  part  to  another  of  a  medium  [4]. 

There  are  several  motivations  for  investigating  the  nresent  eigenvalue 
problem.  First,  the  operator  A^  may  be  used  in  r.ar.y  instances  in  models  of 
Interfaced  media  of  mathematical  physics,  as  nreviously  indicated.  Second, 
the  operator  A^,  at  first  sight,  seems  to  conceal  some  "singular"  behavior, 
using  the  terminology  of  singular  perturbation  of  lumoed  systems  [18]. 


Third,  the  introduction  of  the  sitall  oarameter  t  in  models  of  interfaced 
media  may  be  purely  artificial,  in  order  to  obtain  an  approximation  of  the 
problem  at  hand.  Fourth,  a  formal  Laurent  series  expansion  in  powers  of  c  is 
derived  for  some  elliptic  boundary  value  problems  in  [24,29].  It  is  not 
clear  what  the  relationship  is  between  the  terms  of  this  expansion  and  the 
eigenvalue-eigenvector  pairs  of  the  operator  .  In  the  sequel,  these 
observations  will  be  fully  investigated. 

The  eigenvalue  problem  involving  several  perturbed  operators  has 
been  studied  in  the  literature  [12,13,16,17,31],  and  the  references  therein. 
However,  the  spectral  analysis  of  stiff  operators  (using  the  terminology  of 
[24])  has  not  been  investigated  previously.  This  chapter  presents  a  general 
formulation  using  bilinear  forms  to  avoid  possible  cumbersome  and  comnlex 
boundary  and  interface  conditions.  The  chapter  is  organized  as  follc«-7s. 

The  eigenvalue  problem  formulation  of  a  class  of  stiff  operators  involving 
two  bilinear  forms  is  presented  in  Section  2.2.  In  Section  2.3,  the  con¬ 
vergence  of  the  eigenvalues  and  the  corresponding  eigenvectors  as  t  0 
is  investigated.  Several  examples  are  given  to  illustrate  the  results  obtained 
in  this  section.  In  Section  2.4,  a  generalization  of  the  analysis  of  Section 
2.3  to  p  +  1  (p  >  1)  bilinear  forms,  is  undertaken.  In  Section  2.5,  it  is 
shown,  with  the  aid  of  three  examples,  that  some  of  the  results  derived  in 
Section  2.3  are  applicable  to  a  larger  class  of  eigenvalue  problems.  This 
is  accomplished  by  relaxing  some  of  the  assumptions  made  in  Section  2.2. 

In  Section  2.6,  formal  asymptotic  expansion  in  powers  of  e  of  the  eigen¬ 
values  and  eigenvectors  is  discussed.  In  Section  2.7,  two  numerical  examples 


are  solved.  The  first  example  illustrates  the  properties  of  stiff  operators 
of  Section  2.3.  However,  the  second  example  elucidates  the  properties  of 
stiff  operators  of  Section  2.5.3.  Finally,  in  the  last  section,  some 
concluding  remarks  as  well  as  some  extensions  of  the  forthcoming  analysis 
are  given. 

2.2.  Eigenvalue  Problem  Formulation 

In  this  section,  the  eigenvalue  problems  of  a  class  of  stiff  operators 
is  foimiulated.  Let  V,  H  be  two  given  real  Hilbert  spaces  such  that  V  is 
dense  in  H  and 

Al)  the  injection  of  V  into  H  is  compact. 

Let  V*  denote  the  dual  space  of  V.  After  identifying  H  with  H*,  one  has 
A2)  V  C  h  C  V* . 

Let  i  =  0,1  be  two  forms  on  V  such  that  the  following  assumptions 

hold: 

A3)  is  bilinear,  symmetric  on  V 

A4)  a^(>fi,ip)  is  continuous  on  V,  i.e.,  there  exists  6^  such  that 

1  B.  wev,  V^V 

2 

A5)  where  >  0  and  p^(’)  is  continuous  semi-norm  on  V 

A6)  FgC'^)  +  is  a  norm  equivalent  to 

A7)  a^(v5,^)  =  0  on  C  V,  where  is  an  infinite-dimensional 
subspace  of  V,  i=0,l. 

A8)  If  'i/ 1 — *■  Lq(v)  is  a  continuous  linear  form  on  V,  null  on  Vq,  there 
exists  <P^V  (modulo  V^)  such  that 


(2. 


Let  ,\i)) ,  be  defined  as 

a  (if, 'll)  =  a  +  ea^Cfftp) 

e  u  i 

a-CiPfiij)  =  a^Cififip)  +  aj^(V>,ij/)  •  (2. 

Now  some  important  remarks  clarifying  the  above  assumptions  and 
definitions  are  in  order: 

Remark  2.1: 

It  can  be  easily  seen  from  (A3)  that  a(f,ijj)  are  bilinear, 

symmetric  forms  on  V. 

Remark  2.2; 

From  (A4-A6),  one  concludes  that  a(f,ifi)  are  coercive  and 

bounded  on  V.  In  particular,  for  sufficiently  small  e,  they  satisfy 


(2. 


where  a  (resp.  v)  is  independent  of  e  and  depends  solely  on  (resp. 

Vq,Vj^)  and  the  semi-norms  i  *  0,1. 

Remark  2.3? 

The  bilinear  forms  a(ip,ilj)  define  selfadjoint  operators  [2] 

A^,  A  C  X(V;V*), 

i.e. , 

a^(ip,>p)  *  (  >  ,V  e  V 

*  (  Ai^,)!;  >  ,  V  f  , Ip  ^  V 

where  (  • , •)  denotes  the  duality  cairing  between  V  and  its  dual  v** 

From  the  preceding  remarks,  one  concludes  that  the  snectra  of  A  ,  A  are 


subsets  of  ]R^,  consisting  only  of  the  point  spectrum  [2,15,42]. 

k  k  + 

The  eigenvalue  problem  for  A  is,  then,  to  seek  S  IR  XV  such  that 

G  C  t- 


A 

e 


k  k 


(2.4) 


The  equivalent  variational  formulation  is 


,  k  .  k,  k  , 


y  ^  e 


(2.5) 


Now,  some  well-known  facts  are  summarized  in: 
Proposition  2.1: 


If  (A1-A6)  hold, 

,  €  V  such  that 
e  k=l 


^00  ^ 

then  there  exist  unique  sequences  ^  » 

(2.4) (or,  equivalently,  (2.5))  is  satisfied.  Furthermore 


1)  0  <  Y  <  Y  <  .  . . ,  lim  Y 
e  —  e  —  ,  ,  £ 

k  -*■  +® 


+«> 


2)  is  a  complete  orthonormal  set  in  H,  i.e.,  in  particular, 

k  I  kH 

(x  »X_)  =  5  (Kronecker  delta) 

£  C 

3)  The  multiplicity  of  each  eigenvalue  is  finite 

4)  The  eigenvalues  satisfy  the  following  minimax  formula: 


Y  =  rain  max  a  (x,x)  (2.6) 

W  C  V,  dim  W=k  x  ^  V,lt  xlj^=l 

Proof:  See  [2,15,42] 

P.emark  2.4  : 

It  is  noteworthy  to  mention  that  Prop.  2.1  is  valid  for  any  positive 
value  of  £,  e.g.,  e  =  1.  In  this  case,  the  eigenvalue-eigenvector  pair 


of  A  is  obtained. 


This  section  starts  with  a  series  of  lemmas  which  characterize  the 


various  properties  of  the  spectrum  of  A^.  Then  the  convergence  of  the 
eigenvalues  and  their  corresponding  eigenvectors  as  e  ->  0  is  stated  and 
proved  in  Theorem  2.1.  Some  typical  examples  are  given  at  the  end  of  this 
section,  to  illustrate  the  ideas  advanced  in  the  course  of  the  present 
analysis . 

One  way  to  gather  information  about  the  behavior  of  a  single  eigenvalue 
as  £  -►  0,  is  to  bound  it  from  below  and  from  above  by  known  functions  of  e. 
This  task  is  accomplished  in: 

Lemma  2.1: 

For  sufficiently  small  positive  e,  the  following  estimate  holds: 


k  k  k 
1  Yg.  1  Y 


(2.7) 


for  k  -  1,2,...,  where  {y  are  the  eigenvalues  of  the  operator  A,  i.e., 

they  satisfy 

+  a^(p^,Vj)  =  y'^(p^,V>),V  ^  e  V 
Proof:  For  sufficiently  small  p ,  one  has 


ea(^,i^)  _<  ■<  a(s(>  ,>p)  ,  ¥  ^  f  V  .  (2.8) 

Using  the  mlnimax  characterization  of  eigenvalues  (2.6),  one  readily  deduces 
(2.7)  from  (2.8). 

Now  an  upper  bound  for  the  eigenvector  norm  in  V  is  derived: 

Lemma  2.2: 

Ic  ] 

If  is  any  normalized  eigenvector  of  A  ,  corresponding  to  , i.e., 


/ 


for  k  =  1,2, .. . 


Proof ;  For  each  k,  the  sequence  \  is  bounded  by  Lemma  2.1.  Let  = 


in  (2.5)  to  get 


,  k  k.  k,  k  k. 
a^Cx^.X^)  =  Y^(X^,X^) 


=  Y 


1  Y 


ITot?  one  easily  ^ets  (2.9)  by  using  (2.3).  At  this  point,  the  tools  necessary 
for  finding  the  limits  of  the  eigenvalues  are  available. 

Lemma  2.3: 

The  sequence  '(y^).  is  decomposable  into  t%-70  subsequences  ,  (m'*) 

such  that,  for  each  k. 


lira  =  0  (2.10) 

e  ^  0 

lira  >  0  (2.11) 

e  ^  0  ^ 


k  =  1,2, . . . 

Proof ;  The  following  three-step  contradiction  argument  is  used  to  ascertain 
the  above  lemma. 

1)  Suppose  lim  •  Let  ^  »  where  the 

orthogonality  is  that  of  V.  Take  v  €  and  write  it  as 


V 


Z 

k=l 


f  hv  k 

(v,x^)x^ 


(2.12) 


L’ith  this  choice  of  v,  the  following  inequality  holds: 


(2.13) 


'C  i 


^  ■  y  " 


a^(v,v)  =  aQ(v,v)  +  ea^(v,v)  _>  C 


for  sone  strictly  positive  constant  C,  vjhich  is  independent  of  e. 


Using  (2.12),  a^(v,v)  -  I  Yg(v,x^)^,  which  converges  to  zero  as  e  •>  0, 


contradicting  (2.13). 


k  k 


2)  Suppose  lim  =  Yg  >  0,  k  =  1,2 .  Select  v  €  v  and  write 


e  -»■  0 


it  as  in  (2.12).  It  is  clear  that  the  following  inequality  holds: 


a  (v,v)  =  E  Y^v,xS^  >  C 
k=l  ^  ^ 


(2.14) 


for  some  strictly  positive  constant  C,  which  is  independent  of  e.  However 

^g.(v,v)  =  eaj^(v,v),  which  converges  to  zero  as  e  ->  0,  contradicting  (2.14). 

k 

3)  Suppose  lim  =  0,  for  k  =  l,2,...,i,  without  loss  of  generality, 
e  -*•  0  ' 

Let 

V  =  span{x\x^,* . 

*-  e  e  e 

Select  v  €  v|j  n  v|  and  go  to  step  2,  to  conclude  that  is  infinite 

dimensional,  unless  is  degenerate  (i.e.,  finite-dimensional).  An 

identical  argument  can  be  advanced  to  contradict  the  possibility  that  (2.11) 

is  true  for  k=l,2,...,Z  (£  finite). 

i^j^co  kk®®  k 

Now  decompose  {y  ,x  }i,.i  into  {X  }  if  lim  y  =  0  and  into 

,  othein-7ise. 
e  e  k=l’ 

Some  attention  must  be  focused  on  how  converges  to  zero,  k  =  1,2,... 
Lemma  2.4: 

k  k 

The  sequence  converges  to  zero  with  a  constant  rate  X^ ,  i.e., 

=  1,2,...  . 


•«*  ...W  fc'*  *_*  x'«J!  ‘w*_  ^ 
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Proof :  By  Lemma  2.1,  one  may  assume,  without  loss  of  generality,  that 


,k  -vk  V  .  ,  . 

=  A^e  +  o(e) 


(2.15) 


with  V  €  (0,1].  In  order  to  complete  the  proof,  it  suffices  to  show  that 


Ic  Ic  Ic 

If  ^  is  a  normalized  eigenvector,  i.e.,  ,‘fi  )  =  1,  corresponding 

£  £  £  n  ' 

.  k  , 

to  A  ,  then 
e 


,k  ,  k  k.  ,  ,  k  k. 


(2.16) 


from  which  one  observes  that 


1  O(c^) 


Therefore,  3^(1^^,!^^)  -»•  0  as  e  ->■  0,  which  implies 


(2.17) 


PqCv’^)  ^  0,  as  e  -»■  0 


by  (A5) .  Using  (A6) ,  one  has 


>  C 


(2.18) 


for  some  strictly  positive  constant  C,  which  is  independent  of  e. 
Suppose  that  v  <  1.  Then,  from  (2.15-2.18)  one  concludes  that 


,k  V  ^  . 

A  =  A,e  +  o(e) 
el 


=  ^0^'^e’'^e^  +  o  (e)- 


(2.19) 


From  this,  there  exists  an  element  of  V,  ^  = 


,  k  k  — k,  ,  f 

^1  “  ^0^'^e  +  o  (£  ) 


such  that 


However,  such  a  clain  is  false  because 


— k  — k  1 

e 

— kii 

Since  the  injection  of  V  into  H  is  continuous,  ->■  +“  as  e  0.  In 

conclusion,  there  is  no  element  of  V  such  chat  (2.19)  is  satisfied. 

Remark  2.5: 

k  1 

For  \j  <  1,  the  major  contribution  to  is  supplied  from  (cf.  (2.19)) 
but  the  norm  of  the  contributor  is  concentrated  on  (Cf.  (2.18)),  which  is 
the  paradox. 

Hereafter,  the  focus  will  be  on  the  asymptotic  behavior  of  the  eigen¬ 
vectors.  The  following  lemma  summarizes  the  norm  bounds  of  the  eigenvectors 
Lemma  2.5: 

k  ^  00  k  k 

Let  {A  ,(^  }  ,  ,  be  as  in  the  proof  of  Lemma  2.3,  with  the 

£  £  liC“i  £  £  K”  1 

eigenvectors  normalized  in  H.  Then, for  sufficiently  small  e,  the  follov7ing 
estimates  hold: 


1)  1  (2.20) 

2)  /e  <  C9  (2.21) 

e  V  —  z 

k  =  1,2, .. . 

where  C2  denote  constants  independent  of  t. 

Proof :  Use  Lemmas  2. 2-2, A, 

The  forthcoming  theorem  is  the  main  result  of  this  section.  It  states 
the  convergence  of  the  eigenvalues  and  their  corresponding  eigenvectors 
as  e  ^  0, 
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Theoren  2. 1 

Let  ly  }  be  the  eigenvalues  of  A  and  { y  }  the  corresponding 
normalized  system  of  eigenvectors.  Then,  given  a  sequence  of  e  converging 
to  zero,  {y  ,x  decomposed  into  tv7o  subsequences  {A  ,s^ 

£  £  iiC“J*  £  £  k“i 

k  k 

which  have  the  following  asymptotic  properties,  for  each  k, 

k  k  k 

1)  A^  0  linearly  in  e,  strongly  in  V 

2)  ^  >  0,  ^  weakly  in  H 

e  0  ^ 

k  ®  k 

where  {‘P  and  (i  satisfy 


/  e  Vq  c  V,  vx  e  Vq 


(2.22) 


agC’^^.x)  =  Uq(/,x),  /  e  c  h,  X  e  V 


(2.23) 


Proof ;  Using  the  fact  that  || “  1>  the  estimates  (2.17),  (2.20),  one 

k  k 

concludes  that,  given  a  sequence  of  e  converging  to  zero,  tp  ^  ^  tp  weakly 

in  V  (hence  strongly  in  H  by  compactness).  From  (2.17),  it  results  that 

k  k  k 

<p  S  Vq.  By  Lemma  2.4,  A^  is  asymptotically  equal  to  A^e.  Hence  (2.5) 

degenerates  into  (2.22)  in  the  limit, 
k  k  k 

How  let  w  =  ‘P  -  >p  which  satisfies 

e  e 

,k  _  ,k 

1  ,  k  k^  ,  /  k  k^  e  1^  ,  k  k^  ,  ,  k,  k  k^ 

—  a„(w  ,w  )  +  a,  (w  ,w  )  =  -  (j  ,w^)  +  A  (w  ,w  ) 

e  0  e’  e  1  e’  e  e  e’  e  1  e’  e 


>  a  »w 


^Iv 


The  right-hand  side  of  the  above  equation  converges  to  zero  as  e  0. 


Hence  lw^| 

£  V 


k  k 

0,  indicating  that 'P  ^  -P  strongly  in  V.  It  is  clear  that 


Ic  Ic 

for  normalized  to  1  in  H.  Since  is  0(1),  using  (2.24)  and  the 
ninimax  characterization  of  eigenvalues,  i.e.,  (2.6),  it  results  that 
is  0(1)  and  a^(t|j^,4/^)  is  0(~)  or  equivalently 


=  0(1) 


-  0(^  ) 


Hence, the  estimate  (2.21)  is  tight.  Therefore 


^1  V  ^  as  e  ->■  0 


Note  that  'i/  also  satisfies 
c 


From  (2.25-2.26),  one  deduces  that  a^iip^yX)  is  bounded  as  e  0  and 

a,r>(/'^,X)  is  0(-^).  Since  ||  „  =  1,  taking  formally  the  limit  as  e  ^  0 

1  £  "  £ll  H 

ve 

in  (2.28)  yields: 

an('f'^»X)  =  Un('('^.x)»  ^  X  ^  V 


where  ip  S  (a  subspace  of  H) 


(2 
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Now  consider  the  following  boundary  value  problen 

agC^.x)  +  ea^(v^,x)  =  V  x  ^  V 

Ic 

which  admits  a  unique  solution  v^  €  V  for  positive  values  of  €.  As  e  -»■  0, 
k.  k 

V  -*■  \p  strongly  in  H. 

k  k  k 

Let  w  =  ij)  -  V  which  satisfies 
€  €  € 

aQ(w^,X)  +  ea^(w^,X)  =  "  'i'^^.x) 

+  (p^  -  V  X  e  V  . 

The  left-hand  side  of  this  equation  goes  to  zero  as  e  0,  implying 
('iJg  -  'i'^^.x)  -*■  0  as  £  ^  0. 

Therefore 

-*•  weakly  in  H. 

Remark  2.6; 

The  weak  convergence  in  Theorem  2.1  cannot  be  improved  in  general. 

This  will  be  illustrated  by  Example  2.4. 

Remark  2.7: 

A  careful  examination  of  the  steps  of  the  analysis  undertaken  in  the 

present  section  yields  the  following  observation:  the  weak  limits 
k  0®  k 

>k-i’  ‘'I'  \-i an  orthonormal  system  in  H.  This  remark  is  of  paramount 
importance  in  approximating  the  solution  of  boundary  value  problems  involving 
the  operator  A^. 

Remark  2.8: 

In  the  sequel,  the  fact  that  | "*■  +°°  as  e  -*•  0  is  referred  to  as  the 
resonant  behavior  of 

e 


Now  some  examples  are  given  as  concrete  illustrations  of  the  above 
abstract  results.  Only  operators  of  order  less  than  or  equal  to  four  are 
considered,  due  to  their  frequent  usage  in  modeling  of  physical  processes 
Let  n  =  Qq  LI  U  S  be  a  bounded  set  in  with  boundary  F  =  U  r 

The  manifold  S  denotes  the  interface  between  and  7.^,  as  indicated  in 

Figure  2.1: 


(a)  (b) 

Figure  2.1  a-b.  Examples  of  interfaced  sets. 


Example  2.1:  A  second  order  operator 
Let  H  =  L^(n),  V  =  HpCn) 


aiC#,*')  -  /„  It:  If,  ‘‘='  •  1 

then  (2.4)  becomes: 


J=1  J  J 


where  A  stands  for  the  Laplacian  in  IR 

V  is  the  unit  normal  on  F  or  S,  outward  relative  to 

In  this  example,  (2.22)  becomes 


which  is  a  Dirichlet  eigenvalue  problem  for  the  Laplacian  operator  in 
The  subspace  of  V  is: 

Vq  =  {  X  €  V  :  Xq  =  0,  <= 

r  h-i"” 

The  conclusions  of  Prop.  2.1  are  applicable  in  this  case.  Hence 

2 

is  a  complete  orthonormal  system  in  L  (Q  ) . 
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Equation  (2.23)  becoraes 


Pq  *0  % 


on  Q, 


(2.32) 


which  is  an  eigenvalue  problen  with  mixed  boundary  conditions  for  the 
Laplacian  operator  in  JIq.  Again,  the  conclusions  of  Prop.  2.1  are  applicable 

in  this  instance,  provided  the  interface  S  is  sufficiently  smooth.  Therefore, 

It  2 

{’];  }  form  an  orthonormal  system  in  L  (H  ) .  It  is  noteworthy  to  observe 

0  ^  u 

k  k  2  t 

that  =  0  because  tti  must  be  orthogonal  (in  L  (Jl))  to  ^  ,  1=1,2 . 

The  subspace  of  H  is  then 

*  (X  €  H  :  Xq  €  H^(nQ;rQ),  ^  =  0  on  S,  =  0}  . 


Remark  2.9*. 


In  Example  2.1,  one  can  consider 


(tf  i)  =  E  E  [  a^  dx  i 


where  a,,  satisfies 

jk 


1)  E 

n  n  . 


j=l  k=l  k=l 


Y  r  G  TR^  ,  Z  ^  0. 


The  discussion  therein  remains  unchanged. 
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Equation  (2.23)  becomes 

2  k  k  k 
A  on  Qq 


=  0  on 


(2.35) 


0  on  r. 


k 

^^0  ’  ~  =  0  on  S 


Identical  comments  to  those  of  Example  2.1  can  be  made  here,  provided  some  of 
the  function  spaces  are  changed,  to  reflect  the  increase  in  the  operator  order 
from  two  to  four,  as  seen  in  (2.33). 

For  simplicity  considerations,  a  one-dimensional  version  of  Example  2.1 
is  studied  in  the  sequel.  This  example  will  be  useful  for  illustrating  later 
developments.  It  will  clarify  many  aspects  of  the  eigenvalue  problem  at  hand, 
such  as  nonanalyticity  and  oscillatory  behavior  of  "flattening" 


00 

(and  sometimes  "attenuation")  of  • 

Example  2.3:  (Cf.  Example  2.1) 

Let  fig  =  (a,b)  U  (c,d),  =  (b,c) 


Tg  =  {a,d},r^  =  {d},  S  =  {b,c} 


with  a<b<c<dasin  Figure  2.2. 


0  s  n, 

U  I  1 


Figure  2.2.  Structure  of  the  set  C  ]R. 


^  •/  "k’  0*0  O  *.■■  •.  *.*  •  •.*  /  *.  •-*  *.  -  •-  '•  • 


In  this  instance,  (2.31)  can  be  readily  solved,  to  get 


.k  2  .  ,  x-b 

'P ,  =  — —  sin  kn  — ^ 

1  c-b  c-b 


(2.36) 


The  solution  to  (2.32)  is 
k  '(2k-l)  f' 


fsin(2k-l)  j  ~~  on  (a,b) 


on  (c,d) 


((2k-l)  |)Y 


on  (a,b) 


sin(2k-l)  on  (c,d) 

L  a— c 


=0 


on  (b,c) 


(2.37) 


Remark  2.10: 


Note  that  in  the  case  when  is  not  a  connected  set,  each  subset  of  it 
is  associated  with  a  subgroup  of  eigenvalue-eigenvector  pairs  as  indicated 
in  (2.37). 

Using  Example  2.3,  some  exact  eigenvalue-eigenvector  pairs  are  constructed 
for  some  specific  values  of  e  to  shed  more  light  on  the  asymptotic  behavior 
of  the  spectrum  of  the  operator  .  For  simplicity,  consider  the  following 
example : 

Example  2.4;  (cf.  Example  2.3) 

Let  a=-l,b=Q,  c=d=lso  that  (2.36-2.37)  become: 

=  (kir)^ 
k 

%  -  0 
kl 

s?  =  v2  sin  kTTX 

Mq  =  ((2k-l)  J)^ 

=  cos(2k-l)  y  X 
=  0 

Direct  computation  of  the  eigenvalues  yields  that  must  satisfy: 


Despite  the  transcendental  nature  of  (2.38),  it  is  possible  to  solve  for  y ^ 

for  some  sequence  of  e.  For  example,  for  ,  £  €  D'  , 

2^  (4  f.+l ) 

(u  » }  given  below  is  an  exact  eigenvalue-eigenvector  pair 

(2.39) 

Remark  2.11; 

It  is  instructional  to  observe  the  oscillatory  behavior  of  ^  as 

1  ^ 

il  ->  +».  This  is  consequential  to  the  weak  convergence  of  ^  in 

2  1  2 
L''(fi.)  to  zero,  i.e.,  (ti<  ,  -»■  0  as  2.  V  x  €  L  (n. ) ,  which 

^  L^(n^)  ^ 

can  be  easily  verified.  Furthermore,  it  can  be  easily  seen  from  (2.39) 
that  {i{/  }  are  not  analytic  functions  of  e,  in  the  vicinity  of  e  =  0. 

C 

Therefore,  the  task  of  finding  these  eigenvectors  (in  order  to  solve  boundary 
value  problems  involving  the  operator  A^)  is  nearly  impossible.  The 
alternative  is  to  use  the  weak  limits  of  i.e.,  ^k=l’ 

since  they  possess  some  desirable  properties  (Cf.  Remark  2.7). 

Remark  2.12: 

The  computation  of  ,  using  the  Rayleigh  quotient  of  v 


,  yields 
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1 

M 

e 


£ 


a  ,’4i^  ) 

I  I 

('pI  .pI  ) 

.  £  £ 


fO  ,  TT  .  TT  ,  2,  . 

J-l  (-  2  2’^’  ‘"‘■""j  /o 


2ve, 


sin 


2/e. 


x)  ^dx 


—  +  — 
2  2 


1  /7rv2  ,1  ir  . : 

£ 


2  '2"  ‘2  “£^2/e 


(1)2 

''2'' 


which  shows  clearly  that  the  oscillatory  behavior  effect  is  to  cancel  the 

effect  of  e.  in  the  bilinear  form  in  order  to  contribute  an  amount  of  0(1) 

Z 

to  the  value  of  (Cf.  Theorem  2,1). 

^1  k  k. 

Similarly,  for  e  =  — j  ,  k  =  1,2,...,  {A  ,‘P  }  given  below  is  an  exact 

k  k  k 

eigenvalue-eigenvector  pair 


X  =  (ku)‘ 


e,  0 


/T  sin  /r  kirx  ]  / 

K  K 


e,  1 


=  sin  kTTx 


(2.40) 


Remark  2.13  : 

Note  that  in  (2.40),  depends  on  the  index  of  the  eigenvalue.  Hence, 
one  cannot  let  go  to  zero  and  observe  what  happens  to  this  eigenvalue- 


eigenvector  pair. 


Remark.  2.14: 


The  results  of  Theorem  2.1  show  that 


eO 


0  in  the  present  example. 


This  is  certainly  reflected  in  (2.40)  by  the  presence  of  the  factor 
k  k 

multiplying  ^  _.  This  behavior  of  is  referred  to  as  the  "attenuation" 

e,  u  e 

k  ^ 

of  It  depends  generally  upon  the  choice  of  the  space  V.  See 

Example  2.5. 

Remark  2.15: 

k  k 

The  computation  of  X  using  the  Rayleigh  quotient  of  'f  yields 

^k  % 

,k  f 

— I,  • 

2  k  2 


k  k 

Observe  that  the  relative  contribution  of  a„('r’  t'P  ),  i.e., 

U  e,  e, 
k  k 


,  k  k  . 

(/  )  2 

^k  \  L  (fig) 


,  is  0(£j^).  This  property  seems  to  be  of  a  general  nature. 


See  Example  2.5. 

Remark  2.16; 

k  k 

In  the  proof  of  Lemma  2.4,  it  is  mentioned  that  )  <  0(e),  which 

U  E  e  — 

is  certainly  illustrated  by  the  present  example. 

Remark  2.17 

The  presence  of  /ej^  in  the  argument  of  the  sine  function  in  (2.40) 

indicates  that  "flattening"  occurs  in  as  e  ->•  0  which  is  predicted  by  the 

k  k 

fact  that  ^O(e). 

The  next  example  shows  the  presence  of  flattening  without  attenuation. 


Examp  le  2,5.' 


Let  H  =  V  =  with  2^  =  (-1.0),  =  (0,1),  =  i-: 


=  {!},  S  =  {0}  and  identical  bilinear  forms  to  those  of  Example  2.1 


From  (2.22),  the  limit  of  {s’  },  ,  as  e  0  satisfies 

e  K.=i 


'^0  “  ^0  ' 


dx^ 


,k  k 

\  ^ 


on 


"i  ( 


dV’i 


dx 


KO)  = 


0,  ^J;(1)  =  0 


whose  solution  is 


= 


((2k-l)|-)^ 


0 

.k 


^  S 


cos(2k-l)j  X 


For  E,  =  - j  ,  k  =  1,2, _ is  an  eigenvalue-eigenvector  pair  is 


(2k-l) 


=  ((2k-l)f)2 
k 


v?  =  cos  /e,  (2k-l)—  X 
e,  k  2 

k 


=  cos  (2k-l)-i^  X 
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It  is  noteworthv  to  observe  that  n  in  (2.42)  is  not  multiplied  bv  /e, 

e,  0  k 

k 

which  indicates  no  attenuation.  However,  the  argument  of  the  cosine  function 
does  contain  signaling  flattening. 

Remark  2.18 

Ic  •  Ic 

The  computation  of  X  using  the  Rayleigh  Quotient  of  'f  shows  that 

k  k 


^((2k-l)  4((2k-l) 

1  +  1 
2  2 


Ic  Ic  Ic  Ic  1 

Xote  that  a.(s5  )  =  0(e  )  and  {'P  _)  „  =  —  ,  in  contrast  to  that 

^  ^k  \  k-^  L‘^(nQ)  ^ 

of  Example  2.4.  However,  the  relative  contribution  to  the  value  of  the 

eigenvalue  is  of  the  same  order  of  magnitude. 

In  summary,  flattening  of  is  inherent  to  the  problem,  but  attenuation 

Ic  k 

depends  upon  the  order  of  magnitude  of  which,  in  turn,  depends  upon 

Ic  Ic 

the  choice  of  the  space  V.  If  ^  0(e), then  attenuation  is  present. 


2.4.  General  Results 

In  this  section,  the  results  of  Section  2.3  are  generalized  to  p+1 
(p  >  1)  b’ilinear  forms.  The  generalization  of  Theorem  2.1  is  stated,  bu' 
its  proof  is  omitted,  because  it  follows  similar  steps  to  that  of  Theorem  2.1. 
A  corollary  is  then  presented,  which  considers  positive  as  well  as  negative 
powers  of  e  . 

Suppose  one  is  given  two  Hilbert  spaces  V,  H  as  in  Section  2.2.  Let 
a^iPyV),  i  =  0,1,2,. ..,p,  be  p+1  (p  >  1)  forms  on  V,  with  each  of  them 
satisfying  (A3-A5) . 
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It  is  assumed  that ; 


(A6')  E  p.(i^)  is  a  norm  equivalent  to  j|s5|,, 
i=0 


is  null  on  C  V,  {0} 


(A7')^ 


restricted  to  V^,  is  null  on  V^, 


'a^Cvs,';’)  restricted  to  null  on  /  {0} 

(A8’)  If  liij— *  L.(<j^)  is  a  continuous  linear  form  on  V,  null  on  V  ,  there 
J 

exists  ^  ^j-1  ^ j  ^  such  that 

ajCv’.ii;)  =  L^(tp),  Vii;GV^_^,  j=0 , 1, 2 , .  .  .  ,p  ,  with  V_^  =  V. 

Let  a^(v>,if)  be  defined  as 

P  i 

=  I  e  a  . 

i=0 


Remark  2,19; 

Remarks  2, 1-2. 3  are  applicable  in  this  case. 

The  variational  formulation  of  the  eigenvalue  problem  is  to  seek 
^  IR'*'  ^  V  such  that 

e  £ 

a^(x^,^)  =  ,  Y  e  V  .  (2.43) 

In  this  case,  the  conclusions  of  Proposition  2.1  are  valid.  Moreover,  some 
additional  properties  are  summarized  in  the  following  theorem,  which  is  the 
counterpart  of  Theorem  2.1: 

Theorem  2.2: 

Ic  Ic 

Let  {y  be  the  eigenvalues  and  {y  the  corresponding  normalized 

£  tC~  X  £  iC”  X 

eigenvectors,  as  derived  from  (2.43).  Then  given  a  sequence  of  c  converging 
k  k  *  ** 

to  zero,  . Xj. can  be  decomposed  into  p+1  groups 


Ic  Ic  Ic 

{u  ,'p  »  r=0,l,  .  . .  ,p-l  and  {X  ,'fi  _ ,  with  the  following  asymptotic 

£  j  -4  C  y  X*  tC—  X  C  £  liC“  X 


properties : 


1)  =  X^'c^  +  o  (e^)  ,  ->■  strongly  in  V 


o\  ^  k  r  ,  .  r.  k  ^  ,k  i  i 

2)  r  =  +  o(£  ),  t{i^  r  '^r 

r=0,l, . . . ,p-l 
Ic  Ic 

where  {v?  and  r=0,l,  •  ■  •  ,p-l  satisfy: 


a^(/,x)  =  ^  \  .  vx  e  Vp 


(2.44) 


k, ,  k 


a  (ip  ,x)  =  U  .X),  ti^  SH  CH,  Vx^V  . 
r=0,l, . . . ,p-l. 


□  (2.45) 


Let 

n 

h  =  E  (vJ.iJ^),  n,m  S  (m+n  >  1)  where  b^(i(?,-), 

i=-m 

i  =  -m,-m+-l, . . .  ,n  are  m+n+l  forms  on  V,  satisfying  (A3-A5) ,  (A6’-A7')  then 
one  has: 

Corollary  2.1 

^00  Ic 

Let  {v  }  be  the  eigenvalues  and  {y  the  corresponding  normalized 

£  R®X  £  R'~x 

eigenvectors,  derived  from 


b^(x^,'^)  =  v^(x^,';’)  ,  V  f 


(2.46) 


k  k  ** 

Then,  given  a  sequence  of  g  converging  to  zero,  {v  ,x  can  be 

k  k  k  k 

into  {X  }  and  {p  .  ,'1^  -  -ni,-m+l, .  . .  ,n-l,  vrith  the 

£  £  iC”  X  £  •  <4  £  >  ^  iC*“  X 


decomposed 

following 


properties : 


1)  A  =  Ae  H-ove),!^  strongly  in  '. 

e  n  E 

o\  k  k  r  ^  ,  r.  ,k  ,k  ,  ,  .  „ 

2) u  =ue  +o(e),  'P  -*■  >p  weakly  in  H, 

€  fV  r  e,r  r 

r  =  -in,-m+l ,  .  •  .  ,n*l 
k  ®  iC  ® 

where  {>fi  and  ^  ~  -m»-nH-l, . . .  ,n-l  satisfy: 

b„(Ax)  =  »>^x).  Vx€V___^  (2.47) 

b/4.'^,x)  -  e  C  H,  vx  s 

r  =  -i!i,-m+-l, . . .  ,n-l.  (2.48) 

Proof : 

Let  nt+n  =  p,  a^(ifi  ,ip)  =  i  =  0,1,..., p  and  in 

Theorem  2.2,  to  get  the  desired  results. 

2.5.  Additional  Examples  of  Stiff  Operators 

Several  operators  depending  upon  one  or  more  small  parameters,  possess 

some  of  the  properties  discussed  in  the  previous  sections,  although  they  do 

not  fit  in  the  axiomatization  of  Section  2.2. 

In  this  section,  three  such  examples  are  investigated.  The  differences 

and  similarities  ^^ith  the  operators  of  the  preceding  section  are  highlighted. 

In  Section  2.3,  the  bilinear  form  of  is  assumed  to  be  coercive  on  V, 

2 

i.e.,  a  (ffififi)  >  ae  ,  V  S  V.  In  the  first  example  of  this  section, 

the  bilinear  form  satisfies 

('^,V’)  +  Aili^il^  ^  ae  ,  V  <p  ^  (2.49) 


for  some  positive  constant  A 


2.5.1.  Neumann  eigenvalue  problems 

Let  K  =  V  =  H^(r2)  where  n  =  Ll  U  s  C  ,  with  sufficiently 

^  3  3’^ 

regular  boundary  T  ~  ^1  interface  S.  Let  dx, 

j=l  “i  j  j 

1  =  0,1.  Then  the  eigenvalue  problem  is  expressed  as 


A  ^ 

*XeO  ■ 

k  k 
^e^eO 

on  Qq 

A 

k  k 

Vel 

on 

_ £11 

=  0. 

^  k 
^^el| 

3v 

r  * 

^0 

3v  1 

(2.50) 


^eO  ^el 


0  „  0 

It  can  readily  be  verified  that  this  problem  admits  A  =  0,  'fi  -  constant 

as  an  eigenvalue-eigenvector  pair,  i.e.,  4  <i>  (null  space  of 

the  onerator  A  is  restricted  to  H^(.Q)/??(A  ),  then  its  spectrum  can  be 
e  e 

decomposed  as  in  Section  2.3,  i.e.,  there  are  two  groups  of  eigenvalue- 
eigenvectors  properties  of  Theorem  2.1. 

The  limits  of  the  eigenvectors  as  c  0  satisfy,  respectively. 


Q  -  constant  on  Qq 


,  k  ,k  k 
^^1  = 


(2.51) 


S 


0 
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(2.52) 


k  1 

where  the  constant  in  (2.51)  is  chosen  so  chat  ^  S  H  {9.). 

Remark  2.20; 

00  1  1 
Note  that  €^(^1)  is  not  dense  in  H  (fi).  Therefore,  the  dual  of  H  (P.) 

is  not  a  space  of  distributions.  Hence  (A2)  is  satisfied  in  a  specific 

sense  [30].  The  same  interpretation  applies  to  Example  2.5. 

The  next  example  is  quite  different  from  the  preceding  examples. 

Mathematically,  it  does  not  satisfy  (A5)-(A6).  Physically,  it  arises  in 

Che  field  of  heat  conduction  when  boundary  convection  is  present. 


(2.54) 


Z 

j=l 


9x. 

J 


dx 


^  ^  dr  . 

Problem  2.53  has  no  zero  eigenvalue.  The  limiting  eigenvalue  problem  is 


aQCx^.V*)  =Y^(X^.'^).  V  sff  €  (fi)  (2.55) 

which  has  a  single  zero  eigenvalue  =  0  with  corresponding  eigenvector 

=  constant.  Therefore, there  is  only  one  eigenvalue  of  (2.53)  which 
,  0 

converges  to  y  . 

Now  the  above  discussion  and  the  convergence  of  the  eigenvectors  are 
formalized  in: 

Theorem  2.3: 

Ic  ®  r  k  T 

Let  {y  }  be  the  eigenvalues  of  (2.53)  and  {y  },  _•■  the  corresponding 

€  tC*l  ^ 

normalized  eigenvectors.  Then,  given  a  sequence  of  t  converging  to  zero, 
{y'*,X^},_-,  can  be  decomposed  into  {y^,X^}  and  such  that 

C  C  rC”X  G  C  G  £  iC  X 

1  N  0  r.  0  0  ,  .  „ 

1)  y^  ->  0  ,  Xg.  ^  X  Strongly  in  V 

2)  -*■  >  0 ,  x^  strongly  in  V,  k  =  1,2,... 

k.  k. 

where  {v  ,  satisfy 

*£  G  k=l 

aQ(x^.^)  =  Y^(x^.'^).  V  ^  €  V  (2.56) 


k=0 , 1 , 2 , . .  . 


First  observe  that 


,11  2  _  ii  ji  2 


+  X  Iki!  1  QIq  Ikll  Y»  “o  ^  V  ^  €  V 


(2.57) 


for  some  constant  X  >  0.  For  e  «  1,  the  following  inequality  holds 

agCs^.v?)  +  a^Cv’,'^)  +  X  +  ea^(v’,‘^)  +  X  ||v’||^,  V  €  V 


(2.58) 


T  r  h  k,  <» 

Let  {v  ,p  }  satisfy 

k=0 


aQ(p'^,<^)  +  aj^(p^,(^)  =  v*^(p^,vj)  V  ^  €  V 


Ic  ® 

with  {p  normalized  in  H.  Then  using  the  minimax  characterization  of 

eigenvalues,  one  concludes 


v^  +  X  >  aQ(x^,X£)  +  +  X 

by  (2.57) 


(2.59) 


k  k 

Consequently,  given  a  sequence  of  e  converging  to  zero,  x^  X  weakly  in  V 
(hence  strongly  in  H) ,  where  x  satisfies  (2.56). 


HT  1  ^  k  k 

Now  let  w  =  X  -  X 

e  £ 

ao(w^,vj)  +  caj^iw^.v?)  =  Y^(w^,¥7)  +  (y^  -  y'^Xx^.v’)  +  ea^(x^,*;;),  V  vs  $  V. 


(2.60) 


Let  ^  in  (2.60)  and  use  (2.57)  to  get 


^  i  "5l  X  “o  ■  ”5 


from  which  it  results,  using  (2.60),  that 


11  w^l  ,  -»•  0  as  e  ^  0. 
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Therefore, 

Ic  k 

X  strongly  in  V. 

Now,  a  simpler  version  of  this  example  is  considered  to  illustrate 


the  flattening  of  x  as  e  ->■  0. 

e 

Let  a  =  (-1,1).  r  =  {-1,1}. 


k  k-,“> 


In  this  instance,  direct  computation  of  {y  ,X  },  i  yields 

£  e  k=l 


25,  ,,  .2 

Yg,  =  (liTT) 


X  =  cos  irx  +  —  sin  5,Trx 
e  l-n 


=  ((25,-l)|-)^ 
e  2 


- — -  cos(25.-l)-^  X  +  sin(25,-l)^  x 


5,  =  1,2,... 

Computation  of  the  pair  tYg»Xg}  is  not  trivial.  However,  the  following 
approximation  can  readily  be  found: 


Y°  =  +  O(e^) 

£ 


X^  =  cos  EX  +  sin  £x  +  0(e)  . 


(2.61) 


It  is  worthy  of  observation  that  x^.  flattens  as  e  -»■  0,  to  become  x  =  1 
on  n  in  the  limit. 

Remark  2.21 : 


The  above  analysis  is  unchanged  if 
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(B3)  b^(s5,4')  bilinear,  symmetric  on  H 

(B4)  b^Cvs,!!;)  is  continuous  on  H 

2 

(B5)  b^C'^.'P)  ^  ,  where  >  0  and  q^(‘)  is  continuous 

semi-norm  on  H 

(B6)  q^('r)  ^  »  V  i;)  €  V,  0 

(B7)  qn(v?)  +  q,  is  a  norm  equivalent  to  |lii?|lu‘ 

U  1  n 

In  the  sequel,  the  behavior  as  e  ->•  0  of  the  eigenvalue-eigenvector  pairs 
of  the  following  eigenvalue  problem  is  investigated: 


k  =  1,2, . .  .  ,  V  ¥=  €  V 


(2.62) 


Theorem  2.4: 

k  1  ^  k 

Let  be  the  eigenvalues  of  (2.62)  and  the  corresponding 

normalized  eigenvectors.  Then,  given  a  sequence  of  e  converging  to  zero, 

kki®  kk®  kk®® 

lY  »X  K_i  is  decomposed  into  {X  ,<^  }  and  {y  }  having  the  following 

£  £  iiC~r  C  C  k”i  £  £ 


properties,  for  each  k: 

k  k  k 

1)  0  linearly  in  e,  -+■  strongly  in  V 


(2.63) 


2)  >  0  affinely  in  e,  ^  weakly  in  V 

e  0  e 


(2.64) 


k  ®  k 

where  {■p  and  {i()  satisfy 


/k  ^.,/k  .  ,k,k  .  k^,, 

,x)  +  ,x)  =  ■^j^('<’  ,x)  .  v'  e  v^,  V  X  e 


(2.65) 


aQ('l''^’X)  =  u^(^^,x) . 

k  (  Y.  >) 

s  r  ''  '>1 

bQ(j  ) 

k  -  1,2,... 


€  ,  V  X  e 


(2.66) 


Using  the  simple  arguments  of  Lemma  2.1  and  Lemma  2.3  adapted  to  the  present 

eigenvalue,  one  concludes  that  there  are  two  groups  of  eigenvalue-eigenvector 

Ic  Ic  ^  k  k  ^  k 

pairs  {X  y'P  }  and  {y  ,tp  }  such  that  X  ->  0  linearly  in  e  and 
£  C  C 

k  k 

Pq  >  0.  Adapt  Lemma  2.2  to  the  problem  at  hand  to  get 

ae  I  >^^1  J  +  0e  II  >^^1  J  <  X J  (2.68) 

since  X^  converges  to  zero  linearly,  it  results  that 

lh^|^<C  (2.69) 


for  some  strictly  positive  constant  C,  which  is  independent  of  e,  and,  hence 
given  a  sequence  of  e  converging  to  zero 

k  k 

weakly  in  V. 

From  (2.62),  one  deduces  that 


,  ,  k  k.  2, 

bo('^g..¥’g.)  £  0(e  ) 


since  X^  =  0(c).  Let  ^  ^0  (2.62)  and  take  the  limit  as  e  ^  0 

k  k  k 

(after  dividing  by  e)  to  get  (2.65).  Now  let  -  'f  which  satisfies 

,  k_^k  , 

,  k  k.  ,  k  k.  ,  k  k,  ,  .  e  1^  ,■  k  k.  k,  k  k. 
Oec  le  c  lee  e  ce  lee 


1  all  y  + 


from  which  one  concludes  that 


^  'P  strongly  in  V. 

00  ^ 

For  the  remaining  eigenvalues,  i.e.,  {y  »  one  observes  that  u  ' 

is  a  uniformly  bounded  sequence  due  to  (B4)  and  (2.67).  Hence  using 

Proposition  2.1  one  deduces  that 

ai(/,ij;^)  £0(1),  i  =  0,1  (2.70) 

which  implies  by  (A5-A6)  that 

<  C  (2.71) 

for  some  strictly  positive  constant  C,  which  is  independent  of  e. 

Due  to  (2.67),  the  following  inequalities  hold 

bi(4-^,i(/^)  £0(1),  i  =  0,1  .  (2.72) 

In  fact,  a  better  estimate  can  be  obtained,  by  a  contradiction  argument,  i.e. 

bi(.^^,.i;^)  £  0(e)  .  (2.73) 

Suppose 

=  0(1)  (2.74) 

for  all  bilinear  forms  a^(P,’^),  b^('^’,ti;),  i  =  0,1  satisfying  the  assumptions 
of  the  present  section.  Then  one  may  select  bilinear  forms  (as  in  Example 
2.6)  such  that 


bQ(Xfy)  +  b^(x,x)  =  (v,x) 


Using  (2.62)  and  the  minimax  characterization  of  eigenvalues  to  obtain 


boO^.'-l'g.)  -»■  (>!'  )  =  1,  as  e  ^  0 

k  k 

and  hence  ■*•  0  as  e  -*■  0,  contradicting  (2.74).  Therefore,  (2.74)  is 

valid.  Hence  (2.64)  is  deduced  from  (2.71)  using  the  same  argument  used  in 
the  proof  of  Theorem  2.1.  Mow  expand  (2.62)  formally  in  powers  of  e  (while 
recalling  (2.73)).  The  zeroth  and  first  order  term  yield  (2.66)  for  x  ^  . 

Example  2.6: 

f\  -i 

Let  H  =  L  (fi) ,  V  =  HQ(f2)  where  Q  =  U  U  s  C  ,  with  boundary 
r  =  Tq  U  and  interface  S,  as  in  Figure  2.1. 

Let 

j=i  1  1  j 

b^('^,'j>)  =  v?  'J'  dx  ,  i  =  0,1 
then  (2.62)  is  equivalent  to 


.  k 

k 

k 

k  A 

.  k 

k 

k 

k  .  1 

-  ^  ^Xcl  =  Vtl  ""l  j 

k 

_  n  ,  ^ 

^cO 

r.  '^el 

u 

1 

k 

k 

>^00 

=  ^el 

U 

1, 

* 

on  S 

%0  <1 

dv 

3v  ' 

Equations  (2.65-2.66)  give  respectively 


% 


-  on  j 


t  Ic  • 

=0  =  0 

1  r^  ’  i|s 


^0 


>i'’^  =  0  on 


;  n 

"^olrQ  °  *  "^ils 


k  . 

Uq  =  1 


k  =  1,2, .. . 


Remark  2.22? 


Once  more,  observe  the  flattening  of  {i^  },  -  on  However,  tv  },  , 

z  k=l  0  e  k=l 

oscillate  and  attenuate  concurrently  on  which  ascertain  their 
ameliorated  convergence,  compared  to  those  of  Section  2.3. 

Remark  2.23- 

Remark  2.7  is  also  applicable  in  the  present  eigenvalue  problem. 

2.6.  Formal  Asym.ptotic  Expansions  of  the  Eigenvalues  and  Eigenvectors  oi  A 


It  was  shown  in  Section  2.3  that  {ij;  },  .  are  not  analvtic  functions  of 

£  k=l 

Therefore , they  do  not  have  an  analytic  expansion  in  po'.v'ers  of  c.  Their  node 
of  convergence  (weak  in  K)  has  the  most  profound  impact  in  solving  boundary 
value  problems.  This  shall  be  clearly  demonstrated  in  the  next  chapter. 
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k 

It  was  also  demonstrated  in  Section  2.3  that  {A  },  ,  are  analytic 

e  k=l 

k 

functions  of  e.  Thus,  it  may  be  possible  to  expand  {A  }  as  V7ell  as  their 

C  iC*“  J. 

k 

corresponding  eigenvectors  {<P  }  in  powers  of  £. 

e  tC”  X 

k  k 

In  the  sequel,  a  formal  asymptotic  expansion  of  A^  and  are  derived. 
Then  the  formal  calculations  are  implicitly  given  for  Example  2.1.  Using 
Examples  2, 4-2. 5,  some  of  the  difficulties  associated  with  such  an  expansion 
are  mentioned. 

For  notational  simplicity,  let  the  index  k  of  the  eigenvalues  be  dropped 
and  let 

1  9  9 

(2.75) 

(2.76) 


]  2  2 

A^  =  A  e  +  A  e  +  ... 


0  1 

^  =  ifi  +  e'fl  +  ... 


Substitute  (2. 75)-(2. 76)  into  (2,5)  and  identify  formally  equal  pov/ers  of  t 
to  get: 

Lemma  2.6; 

The  sequences  {A  and  )j^_q  formally  satisfy: 


aQ(¥’  ,X)  =  0,  V  X  ^  V 


a^(V’°,x)  =  ^^(v’^,x)  , ^  ^ 


(2.77) 


^-1 

,  .1  ^  ^  .i+l.,il-i-l  ,  ^  ^ 

»X)  =  2  ^  ('^  ,X)  -  a.,  (-P  ,x),  V  X  G  V 


i=0 

£ 


,  i  .  ~  -i+l.  £-1  .  ..  _ 

a-|  (v'  ,x)  =  2  A  (s3  ,x),  7  X  ^  V 

i=0  ^ 


(2.78) 


£  =  1,2,... 


The  derivation  of  (2.77-2.78)  is  straightforward. 
Now  their  solvability  is  considered  for  Example  2.1: 


Examnle  2.7: 


In  this  case,  (2.77-2.78)  become  respectively 


.0.0 


on 


-A  =  X  on 


,  ^-1 


i+l  £-i-l 


on  n. 


1  -0 

i"o  ’ 

S  3v 

'0  '‘0 
.5.-1 

S 


£ 


£  _  .i+l  £-i  X 

=  L  \  on  n.. 

^  i=0  ^ 


> 


£  n  I  I 

^i|r^  =  '^i|s  ='^o|s 


(2.79) 


(2.80) 


£=1,2,... 


First,  (2.79)  is  identical  to  (2.31)  with  the  obvious  notational  change 
Hence  one  can  solve  uniquely  for  and  It  will  be  shown  that  (2.80) 

is  solvable  recursively.  Let  £  =  1  in  (2.80)  to  get: 
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% 


=  on 


,  1 
’'^0 

1 

'  O 

’  s 

s 

+  x2.; 

on 

’  1  s 

■  "Jls 

(2.81) 


(2.82) 


10  ^  1  /  2 

The  pair  {A  t  is  computed  from  (2.79).  Therefore,  ^  H  (S)  . 

i.  0  V  b 

Hence  (2.81)  is  a  nonhomogeneous  boundary  value  problem  in  H^(.QQ;rQ) 

from  which  one  computes  In  (2.82),  »'^q  known.  Let  zj  G  H^(n^) 


such  that 


’i|r^  ■  ^ils  ■'^0 


and  let  v’ ~  ^  1  ~  ^1  '  Clearly satisfies 


f  ,  -  Is  rl  ,  1  ,  ,  1  1  ,  2  ,0 

(—A  —  A  )'Pj^  =  AZj^  4-  X  +  X 


'P^  =0  =0 

1  ’IS 


(2.83) 


By  the  Fredholm  Alternative  [42],  (2.83)  is  uniquely  solvable  for  in  HQ(f2^) 

0  2 

provided  its  right-hand  side  is  orthogonal  to  in  L  (H^) ,  i.e., 


(AZj^  +  X  z^  +  X  2 


L^(^2^) 


and  consequently. 


(- 

^ 

,  .S5 1  J  9 


Using  Green's  formula,  one  obtains 

1  ^“^l 


,  0  0. 

9 


(2.84 


Using  a  similar  argument  and  recursively  on  £,  one  solves  (2.80)  and  hence 
obtain  all  of  the  terms  in  (2.75-2.76). 

Now  the  yield  of  the  above  calculations  for  Example  2.5  is  presented. 
Example  2.8:  (Cf.  Example  2.4) 

In  this  case,  one  obtains 


=  (kTr)^[e  +  2e^  +  (1  -  |(kTr)^)e^  +  ...] 

e  3 


=  kir(l+x)[e  +  (^^^^  (1+x)^  +  - !))£“  +  ...] 


^  ~  k:rx  +  kTr(l-x)cos  kTx  e  +  ... 


(2.85 


k  =  1,2,... 

Now,  some  remarks  about  this  iterative  process  are  in  order. 


Remark  2.24: 


This  process  yields  a  nonunique  expansion.  For  example,  the  following 
expansion  is  also  given  by  the  same  iterative  process : 


=  (kTT)^[e  -  2e^  -  ■|■(k^r)^e^  +  ...] 


'^eO  ^  kTr(l+x)[e  +  (-  +  l)e^  +  . . .] 


(2.86) 


=  sin  kTTx  +  (1-x)  (kiTcos  k^rx  +  sin  kiTx)e  + 


k  =  1.2,. 


Remark  2.25: 


It  is  not  easy  to  compute  the  terms  beyond  in  the  general  case. 

See  [5,  13,  16,  17,  31,  32]  for  related  topics. 

Remark  2.26: 

The  exact  eigenvalue-eigenvector  pairs  given  in  Example  2.4  for 

particular  values  of  c  suggest  the  following  conjecture:  the  eigenvectors 
Ic  •*  ® 

vP  L  ,  are  not  analytic  functions  of  e  and  therefore  cannot  be  expanded 


as  in  (2.76) - 
Remark  2.27: 

2  1 

In  this  example  as  well  as  in  Example  2.7,  X  and  do  not  depend  upon 
the  choice  of  z^. 

In  the  next  example,  this  observation  is  not  true. 

Example  2.9:  (Cf.  Example  2.5) 

In  an  iterative  process  similar  to  that  of  Example  2.8,  one  obtains 


For  all  computer  runs,  N  is  selected  so  that  an  accurate  plot  is  obtained 

(N  is  indicated  under  each  plot).  The  subroutine  RSG  from  EISPACK  [36] 

with  single  precision  is  used  to  compute  the  eigenvalues  and  eigenvectors 

of  A  . 
e 

Example  2.10:  (Cf.  Example  2.4) 

In  this  example,  the  entries  of  are 


_ h  _ h 

dx  dx 


(2.89 


where 


i,j  —  1,2,...,N— 1 


a(x)  = 


1  if  X  €  (-1,0) 
e  if  X  €  (0,1) 


For  N  even,  K  can  be  written  as 
e 


2  -1 
-1  2  . 


’  •  -1 

-1  1+e  -e 
-e  -2e 


2e  -e 
-E  2e 


(2.9C 


12  1 

The  eigenvalues  X^,  X^,  are  tabulated  in  Table  2.1  for  e  =  0.1, 
0.04,  0.01,  0.001.  The  corresponding  eigenvectors  are  plotted  in  Figures 
(2.3a-2.3d),  (2.4a-2.4d),  (2.5a-2.5d). 


[  I — 1 — j  I  j  I  I  I  I  j  r  i  I  I  j  I  I  I  r 
.8  -8.5  8.8  8.5 

Figure  2.5c.  i(<g  for  c  ■  0.01,  N  ■  100. 
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TABLE  2.1.  for  e  =  0.1,0.04,0.01,0.001 

see 


z 

A^ 

e 

a2 

£ 

1 

u 

c 

0.1 

0.77230 

4.5742 

2.3462 

0.04 

0.36160 

1.3715 

2.4782 

0.01 

0.09650 

0.38677 

2.2153 

0.001 

0.01069 

0.03997 

2.2538 

It  is  clear  that  most  of  the  features  described  in  Section  2.3.  are 

exhibited  in  these  plots.  First,  observe  the  attenuation  of  k  =  1,2 

as  c  -*•  0  on  0  ,  Second,  note  the  oscillatory  behavior  of  as  e  -*■  0  on 
U  z. 

In  Fip.urc  2.5d,  the  cornr-rs  do  not  belong  to  .  .  They  are  inherent  in 

the  Finite  Element  Method  due  to  the  type  of  functions  selected,  i.e., 
i  N-1 

.  Furthermore,  one  may  add  that  the  first  eigenvalue-eigenvector 
pairs  are  computed  more  accurately  than  their  last  conterparts  [17,39]. 
Consequently,  if  the  eigenvalues  of  .A^  are  ordered  ascendingly, 


is  pushed  higher  and  higher  as  e  0  and  hence  computed  less  and  less 
accurately.  Moreover,  since  h  is  fixed,  the  oscillatory  behavior  of 


would  not  be  captured  by  this  approximation,  unless  h  is  made  smaller  and 
hence  increasing  the  order  of  the  matrix  . 


The  following  eigenvalue  problem  is  analyzed  numerically 


d2  k 

,  k  k  k  /  T  ' 

— T-+  X^o  = 


dx 


^2  k 

^  ^rl  k  k  k 

-e  - ^  on  (0.1) 


> 


dx 


> 


XeoW)  ■  *El 


,  k  k 

(0)  =  S-  (0) 


dx 


dx 


(2.91) 


The  finite  dimensional  approximation  of  the  operator  can  be  written  as 


tl  h 

in  (2.88)  where  M  is  given  by  (2.89).  The  entries  of  are: 


(K^).  . 
e  i,j 

_  r  h 

^  ^  Q  dx 

dc^ 

dl: 

+  J 

r 

Q 

0 

i  j  j  1 
dx  +  G 

n  h 

’’h'h 

i.j  =  1 

,2,...  , 

which  can  be 

explicitly 

written 

as 

r 

2 

-1 

ol 

4  1 

0 

-1 

2  -1 

14.1. 

=  f  I 

-1  2  •. 

« 

1 

+ 

h 

/* 

1 

1 

2+2e  c 

(2.92) 
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E!( 


m 


■'■■■: 


0-^ 


Sv 


The  major  difference  between  the  present  example  and  the  preceding  one  is 

k.  k  k.  k. 

that  ->■  Uq  =  1,  and  the  corresponding  eigenvector  weakly  in 


Hq(-1,1)  ,  k  =  1,2,. 

as  £  -*■  0. 


Consequently,  their  oscillatory  behavior  attenuate 


112 

The  eigenvalues  X  ,  u  ,  u  are  tabulated  in  Table  2.2  for 
e  e  e 

£  =  0.1,0.01,0.001,0.0001.  The  corresponding  eigenvectors  are  plotted  in 
Figures  (2.6a-2.6d),  (2.7a-2.7d),  (2.8a-2.8d). 


TABLE  2.2.  for  £  =  0.1,0.01,0.001,0.0001 


£ 

a' 

£ 

1 

u 

£ 

1 

u 

£ 

0.1 

0.6279 

1.6548 

4.555 

0.01 

0.09148 

1.0806 

1.4521 

0.001 

0.01028 

1.0100 

1.0389 

0.0001 

0.001070 

1.0010 

1.0040 

From  the  theoretical  results  of  Section  2.5.3,  it  is  shown  that 


=  (0,simtx)  strongly  in  H^(-l,l) 


=  (simTx,0)  weakly  in  Hq(-1,1) 

2  2  1 
y  1/  =  (sin  2ttx,0)  weakly  in  H  (-1,1) 

£  '  U 

and  this  is  clearly  depicted  in  the  figures  below. 

The  flattening  of  on  is  also  clearly  visible  in  Fi,;<ure.s  2 . 6.i-2  .  hd , 
The  attenuating  oscillatory  behavior  of  on  is  unquestionablv 

documented  in  Figures  2.7a-2.8d. 


S  -■ 


b 


or  £  •  0.01,  N 


Many  of  the  observations  advanced  in  the  preceding  example  are 
valid  in  this  example  as  well. 

Finally,  one  can  easily  estimate  the  orders  of  magnitude  of  the 

Ic  f~ 

attenuation  and  the  oscillation  frequency  of  on  .  They  are  0(/e) 

and  0(— ^) ,  respectively. 

/F 

2.8.  Conclusion 

2.8.1.  Concluding  remarks 

The  eigenvalue  problem  of  stiff  operators  has  been  analyzed  in 
this  chapter,  via  a  general  formulation  using  bilinear  forms,  to  avoid 
the  complexity  of  explicitly  keeping  track  of  the  various  boundary  and 
interface  conditions.  First,  the  intuitive  idea  that  the  eigenvalues  of 
stiff  operators  are  of  different  order  of  magnitude  as  functions  of  the 

II  It 

parameter  e,  is  rigorously  verified.  Second,  many  concealed  features  about 
the  behavior  of  the  eigenvectors  have  been  exposed,  such  as  flattening, 
attenuation,  and  oscillation.  Third,  the  convergence  of  the  eigenvectors 
of  stiff  operators  as  e  0  has  been  investigated.  This  analysis  is  of 
paramount  importance,  because  it  will  yield  insight  into  how  to  approximate 
boundary  value  problems  involving  stiff  operators.  Table  2.3  summarizes 


the  properties  of  some  stiff  operators. 


SUMMARY  OF  PROPERTIES  OF  SOME  STIFF  OPERATORS 
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2.8.2.  Extensions 


The  concepts  used  in  this  chapter  are  very  general.  Consequently, 
several  operators  can  be  constructed  from  the  ones  examined  herein.  They 
can  be  analyzed  using  the  techniques  and  the  concepts  developed  in  this 
chapter.  For  example,  one  may  combine  operators  of  Sections  2. 3-2. 4  with 
those  of  Section  2.5  to  obtain  multiparameter  eigenvalue  problems.  One  may 
also  consider  nonselfadj oint  operators  such  as 


A 

z 


5  3  0,-1 

^  —  a .  .  (x)  r — 

i.j=l  3^i 


u 

r 


-e  -  a!  .  (x) 

i.j=l  5x.  12  OK 


(2.93: 


where  a7.(x)  is  as  in  Remark  2.9  but 
ij 

a^j  (x)  ^  aj^(x),  k  =  0,l. 

The  same  analysis  holds  for  the  selfadjoint  part  of  A  ,  i.e.,  A^ 


The  operator  A^  can  be  written  as 


A  =  A°  +  A^  (2.94: 

e  c  £ 

where 

A°€X(Hj(f:)  ;H"Vc)) 

Ajei:(HQ(  l)  ;L^  )) 


I  i 


given 
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CHAPTER  3 

APPROXIMATION  OF  FORCED  STIFF  SYSTEMS 


3.1.  Introduction 

In  the  preceding  chapter,  the  spectral  decomposition  of  stiff  operators 
was  investigated  and  the  properties  of  their  eigenvalues  and  eigenvectors  as 
functions  of  the  small  parameter  e  analyzed.  One  major  property  which  is 
summarized  below,  is  the  convergence  of  these  eigenvalues  and  their  correspond¬ 


ing  eigenvectors  as  e  0.  It  was  shown  that  the  spectrum  of  A^  is  decomposable 

Ic  k.  -  Ic  k 

into  two  groups  of  eigenvalue-eigenvector  pairs  {  X  ]  ,  (p  )  with  the 

-  ^  k=l  ^  k=l 

following  convergence  mode: 

k  k  k  . 

1)  0  linearly  in  e,  ^  s?  strongly  in  V 

2)  ^  >  0  ,  weaklv  in  H 

e  0  c 

where  V,  H  are  two  given  Hilbert  spaces,  with  V  being  a  dense  subspace  of  H, 
having  a  stronger  topology  than  chat  of  H.  Moreover,  some  of  the  aforementioned 
eigenvectors  are  nonanalytic  functions  of  e.  Consequently,  they  cannot  be 
expanded  in  powers  of  e.  From  these  established  facts,  one  concludes  it  may 


not  be  possible  to  obtain  "strong"  approximations  of  solutions  of  boundary 


value  problems  involving  the  operator  A^. 

In  this  chapter,  the  above  results,  as  well  as  those  discussed  in 
Chapter  2,  are  employed  to  investigate  the  behavior  of  the  solutions  y 
of  the  following  three  abstract  equations: 


=  f 

Xy 

-T^  +  Av  =  f  ,  v(0)=h 

jt  e- t  ■  £ 


(3.1) 

(3.2) 
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3y 

- ^  +  A  y  =  f  ,  y  (0)  =  h  ,  (0)  =  g  (3.3) 

.,^2  £  C  £  Jt 

dt 

The  occurrence  of  (3.1)-(3.3)  is  very  frequent  in  mathematical  models  of 

distributed  physical  processes  such  as  nuclear  reactors,  heat  exchangers, 

chemical  reactors,  fluid  systems,  vibration  systems,  steel  and  glass  processes, 

etc.  Thus,  it  is  important  to  focus  on  them.  A  logical  question  to  ask  is 

the  following:  since  the  eigenvectors  of  A^  are  not  analytic  functions  of  t, 

is  it  possible  to  derive  "weak"  approximations  to  (3.1) -(3. 3)  using  the  weak 

limits  of  the  eigenvectors  of  A  ?  The  answer  is  in  the  affirmative,  but  in 

c 

doing  so  "something"  ought  to  be  lost.  The  major  thrust  of  the  present  chapter 
will  clarify  this  loss  for  each  of  the  boundary  value  problems  (3.1)-(3.3). 

For  elliptic  problems  (i.e.,  (3.1)),  by  appropriately  modifying  the  weak  limits 
of  the  eigenvectors,  one  may  be  able  to  calculate  a  "strong"  asymptotic 
expansion  of  the  solution  of  (3.1).  In  so  doing,  the  formal  results  derived 
in  [24]  are  complemented.  For  evolution  problems  (i.e.,  ( 3 . 2) - ( 3 . 3) ) ,  the 
concept  of  weak  solutions  [23,30]  is  used  to  define  weak  asymptotic- 
approximation  of  the  solutions  of  (3.2)-(3.3).  Hence,  extensions  of  the 
results  in  [23]  are  accomplished. 

This  chapter  is  organized  as  follows.  In  section  3.2,  the  solution 

Ic  -  Ic 

of  (3.1)  is  derived.  The  weak  limits  of  ■  and  ■  ;  ;  are  modified  by 

"  k=l  '  k=l 

adding  to  them  some  appropriately  selected  functions.  The  rationale  behind 
such  modification  is  that  the  modified  limits  become  elements  of  V  and  b.ence 
can  be  used  to  derive  an  asymptotic  expansion  of  the  solution  of  (3.1). 


M  T 


In  Section  3.3,  the  convergence  of  the  solution  of  (3.2)  is  investigated  and  a 
weak  asymptotic  approximation  is  constructed  for  it.  In  Section  3.4,  an 
analysis  similar  to  that  of  Section  3.3  for  a  class  of  hyperbolic  problems 
(i.e.,  (3.3))  is  undertaken.  In  Section  3.5,  some  concluding  remarks  and 
extensions  are  presented. 


3.2.  Elliptic  Boundary  Value  Problems 

As  in  Section  2.2,  suppose  two  Hilbert  spaces  V  and  H  are  given. 

The  same  notation  and  the  same  assumptions  are  kept. 

First,  the  modifications  of  the  weak  limits  of  the  eigenvectors  are 
considered.  Then,  using  these  modified  limits,  an  asymptotic  expansion  of  the 
solution  y^  of  the  following  boundary  value  problem  is  constructed 


A  ,_=  f  ,  fSH 


(3.4) 


or  equivalently 


a^(y^.,'F)  +  £:a^(y  ,v')  =  (f,>;)  ,  y  _ev  ,  Vv^SV 


(3.5) 


From  the  preceding  chapter,  the  weak  limit  of  ■■'j  ^  in  H  and  the  weak 

,  .  .  k  ,,  .  _  .  , 

limit  of  in  V  sati&ry,  respectively: 


/  k  ^  k.  k  ,  k_„  ,, 


(3.6) 


,  v.sv^ 


(3.7) 


k=l,2,3, . . . 


-pec- ivc  ly ,  V^^)  is  a  subspace  of  H  (respectively. 
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Since  ■.  ,  one  may  add  to  it  a  function 


-  rk  ,  _2  ,k 

"0  -1  ^  ^2  •  • 


(3.8) 


where 


H  such  that  G  V 


(3.9) 


V  ,  i=i,2,... 


It  was  noted  in  the  preceding  chapter  that  -  v  i  and  '  s?  )  form  an 

k=l  ,  k=l 

orthonormal  system  in  H.  Therefore,  one  requires  that  satisfy 


aQ('/  +  -^.X)  +  =  Uq(-/,x),  VxGV 


(3.10) 


substitute  (3.8)  into  (3.10)  and  identify  formally  equal  powers  of  ,  to  get: 


a^d^.x)  =  0  ,  V>:€V 


,(-:^,x)  +  a^(lj_^,x)  =  0  ,  VxGV 


(3.11) 


=1,2,.. 


from  which  one  finds  chat 


'■  1=0,  k=l 


a^(:Q,x)  =  0  ,  €  V  ,  VxGVq 


(3.12) 


,  V,GV 


a^(C,  <)  =  0  , 


4=1,2,.  .  . 


(3.13) 


k=l,2,... 


Remark  3.1:  It  is  worthy  of  mention  that  the  iterative  process  described  by 


( 3 . 11) -(3 . 13)  appears  to  average  the  oscillatory  behavior  of  v 
Similarly,  one  adds  to 


_,k  ,k  2,k  . 

3  =  e  o,  +  e  +  .  .  . 

e  1  2 


(3.14) 


where 


3'"  €  V  ,  1=1,2, ■ 


The  function  ■■  is  chosen  such  that 


+  £a^(^^  +  6^,.)  =  £>-^(^^,X)  ,  Vx^V 


Ic  Ic  ® 

from  which  one  concludes  (using  the  fact  that  P  that  (3  } 

^  ^  ?,=l,k=l 

3^(6^,  X)  =  -  a^(^‘',:<)  ,  Vx€V 


a^(-J,X)  =  0  ,  VySV^ 


satisfy 


,.k  , 


‘1 


(•:-,<)  = 


-  ,  v<ev 


0  ,  v,ev 


i=2,3,. . . 
k=l,2,. . . 


(3.15) 


Remark  3.2:  The  zeroth  term  in  (3.14)  is  zero  because  p  ^V. 


6 


Remark  3.3:  The  iterative  process  presented  above  appears  to  average  the 

Ic 

flattening  (and  attenuation)  of  } 

^  k=l 

Before  proceeding  further,  consider  the  following  two  examples 


examine  what  (3.11)-(3. 13) ,  (3.15)  yield: 


Example  3.1:  (Example  2.1  continued)  Let  H  =  L  (fi) ,  V  = 


r  9  J’p 


Hj(n) 


.  ,  i  ox.  9x. 

J  =  1  3  J 


then  ( 3 . 11) -( 3 . 13) ,  (3.15)  become 


Cqo  =  0  on 


"01 

=  0 

on 

k  1 

=  0 

.k 

01k 

’  "01 

""-^0  =  °  "'O 


.  -k 
9v 


. ,  =  0  on  n. 
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The  above  equations  represent  boundary  value  problems  for  the 
biharmonic  operator  in  each  region  i=0,l. 

Remark  3.4:  The  computation  of  these  modifications  is  unquestionably  easy, 
because  the  dependence  on  the  parameter  e  is  eliminated. 

Now  attention  is  focused  into  obtaining  an  asymptotic  expansion  of  the 
solution  of  the  boundary  value  problem  (3.1)  using  the  modified  weak  limits  of 
the  eigenvectors  of  A^. 

It  is  noteworthy  to  indicate  that  (3.4)  (or  equivalently  (3.5))  has 
been  solved  in  a  more  general  context  than  here  in  [24],  Therefore,  only  the 
details  pertaining  to  the  present  approach  are  given. 

In  the  next  theorem,  the  usage  of  the  modified  limits  is  shown. 

Theorem  3.1:  For  sufficiently  small  e,  the  solution  of  (3.4)  is  given  by 

"  c*"  k  k  “  d*'  k  k 

y,  =  2  \  (/  +  ?,)  +  (^  +  ^,)  (3.16) 

k=l  Uq  “  k=l  e'-j_ 

where 


(3.18) 
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k  ^  Ic  X 

Proof :  The  eigenvectors  ^  ^  ^k~l  ^  complete  orthonormal 

Ic  k  k  k 

basis  of  H  and  -  +  S  V,  €  V  by  construction.  Now  it  is  a  matter 

of  verification  that  (3.16)  is  the  unique  asymptotic  expansion  of  y_. 

k  k 

Remark  3.5:  Truncate  C  and  6  to  the  Pth  term  and  denote  these  truncated 
-  £  £ 

k  n  k  D  P 

series,  respectively,  by  C  and  0  Define  e^  by 

£  £  £ 

P  P 

=  V  - 

£  '  £  £ 


where 


cc  Ic 

,P  -  c 


d _  ^  k  k-,p,, 


y.  =  2  ^  (v  +  c/n  +  r  (VJ  + 


k=l 


k=l  £>.: 


Then  it  can  be  shown  ([24],  p.  13)  that 


lle^^ll  <  C  £p 
£  V  - 

where  C  is  a  constant  independent  of  e. 

k  k 

Remark  3.6:  Note,  that  if  the  forcing  f  equals  >p  (respectively,  v  )  then  y^ 

1  k  k  1  k  k 

becomes  — —  +  6  )  (respectively,  — ^  (+  +  )).  This  clarifies  Remarks  3.1, 

.  fC  c  K  £ 

-  P  n 


Example  3.3:  (Examples  3. 1-3. 2)  In  this  case,  (3. 16)-(3. 18)  become 


"  ,,k  .  -k  ^  :  d''  k 

k  *''0  'eO^  ,k  eO 

k-l  k=l  .;.j 


k  k 

V  <-  -k  .  ^  d  ..k  k 

k=l  k=l  £•  ^ 

0  1 


Ik’" 

Note,  '  ,  ,  are  renormal  i-ced  in  H. 

'  Lr  =  1 


where 


k  /p  k, 

c  =  (f.'j;  ) 


2  2 
L  (a)  ^  L^(2q) 


L  (^2)  ^  ^  L^(a^) 


k  k 

are  computed  in  Example  3.1  or  3.2, 


Remark  3.7:  Note  that  lim  II  y  ^  is  finite  but  not  lim  tly^-,11  ^ 

-0  L^Qq)  £-0  L^(2^) 


Remark  3.8:  By  examining  (3.16),  one  concludes  that  the  group  of  eigenvalues 

k 

that  accumulate  at  zero  as  e  ->■  0,  i.e.,  {A  },  ,  causes 

t  k=l 


lim  II  y  .  II  „  -).  +  CO 

c^O 


Remark  3.9:  The  eigenvalue  problems  considered  in  the  preceding  chapter  can  be 

interpreted  as  the  boundary  value  problem  (3.4)  with  f  depending  nonanalytically 

k  k 

on  e,  i.e.,  f  =  f(x,£)  =  However,  it  can  readily  be  seen  that  the  process 

by  which  (3.4)  is  solved,  is  not  applicable  in  this  case. 


Remark  3.10:  If  f  were  analytic  in  s  and  hence  expandable  as 


0  1  2  2 
f  =  f  +  ef^  +  t  f  + 


then  using  the  above  procedure,  one  can  solve 


A^_y^  =  f"^  ,  r=0,l,2. 


to  get  an  approximation  of  the  solution  of 


>> 
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3.1. 

{1}. 


1+-: 


(1+x)  J 


eO 


=  V 


lc=l  (.^) 

:  2 1^  ^ 
k^l  e(>^)'/2  1+c  1 


(  1+x) 


_  ,  2(-l)^'^  1  ,,  , 

^£l  ,  ",  .  k,3/2  1+r 


k=l  (Uq) 


,  2 

k:i  £(aJ)^/2 


R  .-  /k 

[sin  A  X  +  - - A.  (1-x) 

1  l  +  £  1 


(3.19) 


where  no  simplification  is  made  to  allow  the  origin  of  each  term  to  be  identified 
The  exact  solution  can  be  computed  directly  and  is  given  by 


X 


_+  1  iz£ 

2  2  1+e 


X  + 


1 

l+£ 


-;1 


x"  1  1-t  .  1 

2t  2  e(l+£)  ^  1+t 


> 


I 


(3.20) 


It  can  be  easily  shown  that  (3.19)  is  identical  to  (3.20). 

Remark  3.11:  The  generalization  of  the  techniques  presented  in  this  section  to 
include  p+1  (p>l)  bilinear  forms  is  straightforward  (Cf.  Section  2.4). 

Remark  3.12:  Similarly,  the  techniques  of  this  section  can  be  applied  to 
boundary  value  problems  involving  some  of  the  operators  discussed  in  section  2.5. 


3.3.  Parabolic  Boundary  Value  Problems 

In  this  section,  an  evolution  problem  of  parabolic  type  is  considered. 
Hence  let  the  variable  t  denote  time.  It  is  assumed  that  t  €  (o,T),  T  ■  and 
that  all  the  assumptions  made  in  Section  2.2  hold.  Let  L“(0,T,V),  L‘’(n,T;H), 
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L  (0,T;V*)  denote  the  Hilbert  spaces  of  Lebesgue  square  integrable  functions 
with  values  in  V,  H,  V*,  respectively  (Cf .  Appendix).  Let  prime  denote  the 
distributional  derivative  with  respect  to  time  [30]  .  In  the  sequel,  the 
following  parabolic  boundary  value  problem  is  analyzed 


p  fiD-fli42  364 
1  UNCLRSSIFIE 

RNRLVSIS  RND  CONTROL  OF  R  CLRSS 
DISTRIBUTED  SVSTEMSCU)  ILLINOIS 
RND  CONTROL  LRB  H  SRLHI  JUL  83 

) 

OF  STIFF  LINEAR 

UNIV  AT  URBANA  DECI5I0^ 
DC-61  Neeei4-79-C-e424 
F/G  12/1 

2/2  1 

1 

1  I 

_ J 

^■1 

■HI 

■■■ 

_  ^ - - 

'.he  preceding  discussion  is  summarized  in 


Theorem  3.2:  Let  denote  the  solution  of  (3. 21)-(3. 23) .  Given  a  sequence 
of  e  converging  to  zero, 

-►  y  weakly  in  L^(0,T;H)  (3.34 

where  y  is  given  by  (3.31).  Moreover, 

By  II  i  C  (3.35 

^  L^(0,T;H) 

C  is  a  constant  independent  of  e. 

Proof ;  Use  the  eigenvalue-eigenvector  pairs  of  the  operator  and  their 
properties, as  described  by  Theorem  2.1, to  obtain  (3.34).  The  estimate  (3.35) 
is  then  readily  derived  by  employing  (2.21). 

Remark  3.13:  The  convergence  of  y^  in  Theorem  3.2  cannot  be  improved  (cf. 

Jq  00 

Remark  2.7).  Also,  it  should  be  emphasized  that  }  in  (3.31)  are  not 

k=l 

renormalized  in  H. 

3.3.2.  Asymptotic  approximation  of  y^ 

It  should  be  noted  that  the  method  by  which  stiff  elliptic  boundary 
value  problems  were  solved  in  section  3.2  does  not  yield  an  iterative  process 
for  evolution  problems  in  general.  Therefore,  one  would  be  content  to  obtain  a 
"weak"  approximation  of  the  solution  of  (3. 21) -(3. 23)  using  only  the  weak  limits 
of  the  eigenvectors  of  A  . 


In  the  sequel,  only  a  more  general  version  of  Example  2.1  is  considered, 
i.e.,  the  operator  is  written  as 


A  0 

A  =  °  ,0<e<^l 

^  0  eA^ 


where 

"  "  k 

A,  =  E  E  a  .  (x)  ^  ,  k=0,l 

K  .  ,  ,  ,  OX .  11  OX. 

1=1  J  =  1  J  J 

k 

a^^ (x)  satisfy  the  conditions  of  Remark  2.9.  The  chief  reason  for  this  digression, 
is  that  one  can  explicitly  specify  the  regularity  conditions  of  the  functions 
involved  in  the  construction  of  the  approximation  of  y^.  However,  the  concepts 
involved  herein  are  equally  applicable  to  the  general  case. 

First,  the  following  notation  is  adopted  in  the  sequel. 

Q.  =  X(0,T)  ,  i=0,l 

E.  =  r.  X(0,T)  ,  i=0,l 
11 

R  =  S  X(0,T) 

0  0  0 

Now,  let  the  zeroth  order  approximation  be  denoted  by  y  =  ^^0’^  1^ 
and  defined  as  follows.  Let  e  -»■  0  formally  in  (3.21)  and  retain  the  following 


ir  *o>'o  *  'o ‘'o 


(3.36) 


yo  =  0  on  Eg 


( 3. 37) 


0  on  R 


yo(0)  =  Hq  on 


The  solution  of  of  (3. 36)-(3 . 38)  is  regular.  Actually, 

y°  €  L^(0,T;H^(np;rg))  C  L^(Qq).  Hence  y°  e  (S) )  C  l^(R) 

R 

Let  the  second  part  of  satisfy 


X  A  0  f 

ir  "^^el  =  ^1  ^1 


y^l  =  0  on 


0  0  D 

y,l  =  Vq  on  R 


y^l(O)  =  on  Hq  .  (3.4 

Problem  (3. 39)-(3.41)  is  a  nonhomogeneous  boundary  value  problem.  Consequently 
y^^  has  meaning  in  a  weak  sense  using  transposition  [30]  . 

Since  the  zeroth  order  approximation  y^  is  weak,  one  rewrites  (3.36)- 
(3.53)  in  the  proper  form,  using  transposition  as  follows.  Let 


^0  ^0  R 


^  H  ’  (Q^)  ,  =0,  Xj^ 


0  ,  X;l(T) 


Now  consider  the  following  isomorphisms 


■^0  ^^0 

Xq  H-  3^-  +AqXo  from  to  L‘(Qq) 


(3.4 


from  to  L  (Q^)  . 


(3.43) 


By  transposition,  one  concludes  that  Xq  *  M  (Xq).  being  a  continuous  linear  form 

2  1 

on  (endowed  with  the  topology  induced  by  H  ’  (Q^)),  there  exists  a  unique 
Yq  e  L^(Qq)  such  that 


q/o<-  *oXo)  % 


do  =  ^  (Xn)  .  Vx^6$, 


(3.44) 


»  0  I  0  . 

M  1 — *■  yQ  is  a  continuous  linear  mapping  of 


Remark  3.14). 


(3.45) 


Similarly, 


!”*■  being  a  continuous  linear  form  on  (endowed  with  the 

2  1  0  2 
topology  of  H  ’  (Qj^)),  there  exists  a  unique  y^^  €  L  (Q^)  such  that 


I — >■  y^j^  is  a  continuous  linear  mapping  of  1'*  | — ►  L^(Qj^) 


(3.46) 


(3.47) 


00 

Remark  3.14:  Since  (Q^)  cf.  the  dual  of  is  not  a  space  of  distributions. 

Therefore,  the  introduction  of  ©_  is  required  to  interpret  duality  [30]  . 

‘0 


■  O  « .•  'V-  •  •  ■ 


Select  M  (Xq)  and  M  (x^)  as 


where 


M  (Xq)  -  /  I 

Qn 


mJ:(x,)  =  /  fiXidPo  +  /  h  Xi(x,0)di2 
^  ^  Q- 


fi  e  L^(Q^)  ,  h.  6  L^(Q.)  ,  i=0,l 


(3.48 


(3.49 


(3.50 


Theorem  3.3:  Given  f  ,  h  as  in  (3.50),  M^,  as  in  (3.48)-(3.49) ,  there  exists 
a  unique  ^0  ^  ^ei  ^  such  that  (3.44),  (3.46) 

are  satisfied. 


Proof ;  It  is  easy  to  verify  that  the  solution  of  (3.44),  (3.46)  is  identical  to 
that  of  (3.36)~(3.41) .  Let  Xq  ^  Cg  (Qg)  in  (3.60)  to  obtain 

*  ''  fo^o'^^O 
^0 

and  hence 

(  0  ^^0  r  “ 

Q  ^0^"  It”  ^O^O^^^O  *  ^  ^0^0^^%  ’  '‘^O  ^ 

0  Qg 

which  yields  (3.44).  Now  multiply  (3.36)  by  Xn^'^n  and  apply  Green's  formula  to 


dy^  3x 

[  ^oV%  “  ■  i  yoXo^^’^^'^S  ^  I  ^0^'  JT  ^  • 

%  %  ^  ^0  % 

Hence  by  comparing  it  with  (3.44),  (3.48),  one  can  obtain  (3.37),  (3.38). 


Likewise,  one  can  easily  verify  that  the  solutions  of  (3 . 39) -(3 . 41) 


and  (3.46)  are  identical.  Now  an  error  estimate  between  the  exact  solution  of 
(3. 21)-(3 . 23)  and  y^  is  derived. 

Theorem  3.4;  Let  y^  be  the  solution  of  (3 . 21) -(3 . 23)  and  y^  the  solution  of 
(3.44),  (3.46).  Then,  for  sufficiently  small  e,  one  has 

II  y.  -  y°ll  2  ic  (3. 

L^(Q) 

where 

Q  =  Qq  X  Qi 

C  is  a  positive  constant,  independent  of  e. 


Proof :  First,  rewrite  (3.21)-(3.23)  in  the  weak  form,  i.e.. 


^0 


a 


0 


"0  0 


I  ^  =  I  /  h^Xi(x,0)dfiQ  -  t/y^^  ^  dR,  Vx^G^^.  ( 


r> 

“1 


Subtract  (3.44)  from  (3.52)  and  (3.46)  from  (3.53)  to  obtain 


9X, 


(3. 


where  the  following  interface  condition  is  used 


Since  -*■  y  weakly  in  L  (0,T;L  (Q))(-  L  (Q)) 


L^(0,T;L‘-(n)) 


since  y°  €  (0,T;H^(S^q; Fg) )  C  l^(Qq) 

then  y®  e  L^(0,T;H^^^(S))  C  L^(0,T;H~^^^(S) )  .  (3.5 

R 

From  (3.58),  (3.59),  one  deduces  that  ^  continuous  linear  form 

on  0^  ,  i=0,l  (3.6 


0  I  0  *  I  2 

y^Q  -  y^  is  a  continuous  linear  mapping  of  ^  ^^0^ 


iw 


where  C^,  are  some  positive  constants,  indepenaent  of  £.  Hence, 


»y  - 

L^(Q) 


=  lly  „  -  y^ll  , 
■^eO  ^0  ,2 


L  (Qq) 


2 


<  C  e 


Remark  3.15:  One  ought  to  mention  that  (3. 51) holds  for  t  small,  but  strictly 

positive,  because  of  the  heavy  reliance  upon  the  fact  that  the  solution  of 

(3.64)  belongs  to  It  is  easy  to  see  that  if  one  sets  formally  e  =  0  in  the 

2  2 

first  equation  of  (3.64),  its  right-hand  side  would  be  in  L  (0,T;  L  (J^^))  but 

X  for  £  =  0  would  not  be  in  . 

Cl  1 

Recall,  that  it  was  shown  in  Section  2.3  that  the  eigenvalue-eigenvector 

k  k  k  k 

pairs  of  A  ,  i.e.,  {y  ,x  }  are  decomposable  into  two  groups  {X  }  and 

,  ,  ’  ^  ^  k=l  ^  ^  k=l 

{p  } 

^  ^  k=l 

1  k 

The  normalized  weak  limits  (in  H_(f2))  of  satisfy 

■  k=l 


V^O  =  0  on  fig 
.  k  ,k  k 

on 


=  0,  =  0 


(3.65) 


k*l , 2 , . . . 


The  normalized  weak  limits  (in  L  (fl))  of  {<p  }  obey 
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k  ^  k  k 
O'O  ''O  '0  '“0 


"Jq  =  0  on 

-ol  -0. -^1  -0 
°  I  A  I  C 

'  ‘0  0 


J 


k=l,2,... 

Then  the  solution  of  (3.44),  (3.46)  may  be  represented  uniquely  by 


(3.66) 


0  „  k,  ,  k 

’'o  •  "  "^>'-0 

k=l 


00 

^  €  L^O.T),  E  /  |c‘'(t)|^dt  < 
k=l  0 


(3.67) 


0  _  ,k.  ,  k 

y  =  E  d_(t)>^ 

k*l  ■-  ^ 


t,  ")  00  X 

d_  e  L‘-(0,T),  E  /  |d^(t)l^dt  <  00 
'■  k=l  0 


(3.68) 


k  k 

In  order  to  detemine  c  (t)  ,  d  (t)  ,  let 

e 


XQ(x,t)  =  6(t)<iQ(x),  6  €  C^([0,T]),  9(T)  =  0  (  so  that  Xq^^'q) 


:<j^(x,t)  =  v(t)<^^(x),  V  €  C^([0,T]),  ;(T)  =  0  (so  that 
in  (3.44),  (3.46)  to  get 


/  c‘^(  -  +  p^6)dt  =  /  (fo.Vo)-^dt  -  (hQ,l.J)-(0) 


(3.69) 


i 
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T  T 

/  d^i-  +  X^ev)dt  =  /  (f  .s?hvdt 
^  £  dt  1  0  1  1 


(h^.V’pvCO) 


^  CO 

e  ja  C^f  'v^  -;^dS)vdt,  (3.70) 
0  e=l  S  ''^A^ 


which  are  equivalent  to 


dc  ,  k  k  k. 

J^  +  UqC  =  (fg,.4^Q) 


(3.71) 


c^(0)  =  (hQ,<j/Q) 


e  ,  ,k  ,k  k, 

—  +  X^ed^  =  (f^.>^P 


-  e 


CO  5^ 

^  C  j  Ip^  -r -  dS 

e-i  s  ° 


d^<0)  -  (hj,»J) 


(3.72) 


Remark  3.16;  Using  the  terminology  of  singular  perturbation,  one  concludes  that 
(3.67)  represents  the  "fast"  subsystem,  while  (3.68)  represents  the  "slow" 
subsystem  driven  by  the  fast  one. 


Remark  3.17:  Note  that  v  , 

-  -  el 


=  y 


0 


has  meaning  in  an  "average" 


sense.  Hence,  there  is  no  contradiction  in  representing  y^^  on  S  with  the  aid 


of  which  are  null  on  S  [23] 


Remark  3.18:  If  f  =  (0,f,)  and  h  =  (0,h,)  then  y?  =  0  and  hence  v*^  €  L^(0,T; 

i  XU  c 


Hf.(n)).  In  conjunction  with  this  remark,  one  adds  that  if  is  the  empty 


- - -  - - ,  w..-.-  .....  -.Q 

set,  then  problem  (3.21)-(3.23)  degenerates  to  a  problem  analyzed  in  [24]. 


Remark  3.19:  At  this  point,  it  may  not  be  clear  how  better  convergence  of 


y^  -  y^  (as  evidenced  by  the  estimate  (3.51))  is  achieved.  Uae  to  the  weak 


convergence  of  }  in  H,  by  taking  the  limit  as  c  0  in  (3 . 21 ) - (3 . 13) , 

^  k=l 

2 

something  of  0(1)  (in  L“(0,T;H))  is  lost  from  y_.  However,  by  renormalizing 

k  k  ^ 

the  weak  limits  of  }  ,  i.e.,  {ijj  }  ,  this  deficiency  is  corrected. 

^  k=l  k=l 

Intuitively  speaking,  this  action  amounts  to  saying:  since  converges  to 

2  k 

zero  weakly  in  L  (H^) ,  keep  all  of  the  energy  associated  with  'ij ^  in  region  2^. 

It  is  noteworthy  to  report  that  this  phenomenon  does  not  occur  with  , 

'  k=l 

because  they  converge  weakly  in  V  (hence,  strongly  in  H) .  Therefore,  their 

k 

weak  limits  are  automatically  normalized.  However,  a  -*  0  linearly  as  ->  0. 

k  k 

Hence,  a  must  be  replaced  by  its  asymptotic  equivalent,  i.e.,  \  e,  and  not 
^  J. 

its  limit,  to  get  a  better  estimate.  At  present,  consider  the  following  examplt 


Example  3.6:  (cf.  Example  3.4).  Let  f^  =  1,  f,  =  0,  h  =  0.  Then  (3 . 21) - (3 . 23) 


become 


3y  3  y  n 

EQ 


ax 


■|—  =  1  on  (-1,0)X(0,T) 


9 

3“y 


El 


ygQ(0,t)  =  y^^(0,t) 


■^--0  ^>^.1 

- -(0,t)  =  E  — ^(0,t) 


3X 


?x 


at  “  ^  — F  =  0  (0,1)X(0,T)  j 


y_Q(-l,t)  =  y^^(l,t)  =  0 


(3.73) 


f 


y  (x.O)  =  y  (x,0)  -  0 


I 
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Recall  that 


Uq  =  ((2k-l) 


Vq  =  /2  cos  /UqX 


/k 

rU  -J 


xj  -  (kit)^ 


XF» 


sin  v'A^x 


The  solutions  to  (3 . 71)-(3 . 72)  are  easily  computed  and  are  given  by 


e\t)  = 

,  k.3/2 

(Pq) 


k-1 


(1  -  e  ^0*^) 


d^(t)  =  I  2^(-l) 


S,-l 


a=l  ,1.3/2  ,.K1/2 

(Pq) 


(1  - 


i  ,  k 
UQ-^ie 


+  e  s 


2v^(-l) 


l=l 


.-1 

e  0 


Z  ,  k 

V*1'= 


Consequently,  y^  is  written  as 


^0  =  ^ 


k=l  (Pq) 


00  I  CO 


^  r 

^ti  =  ^ 


A(-l) 


£-1 


k=l|£=l  (Pq)^^^ 


(1  - 


'0 


Z  ,k 
Po-\e 


.  k 

^-Aist, 


sin  /A^x 


+  E 
k=l 


e  i 


”  £-1 

,  4(-l)  1  -Pot 

n_i  /,  ^3/2  ^  £  .k 

£-1  (Pq)  Pq-X^c 


£ 

e  0^ 


■  A^ 

sin  /Xj^x 


(3.74) 


(3.75) 


I 


I 


•Ii 

id 

-A, 


.•« 


I* 


M 


I 


I 
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Remark  3.20:  Observe  that  contains  exponential  functions,  decaying  with  rates 
u».  Therefore,  y^  represents  the  fast  subsystem.  In  contrast,  y^  contains  one 
slow  component  and  a  fast  one. 

0  2 

Remark  3.21:  Note  that  y  ,  0  in  L  (Q, )  in  this  case,  which  makes  sense.  See 

numerical  analysis  in  Chapter  5.  ° 

One  may  proceed  further  along  the  same  lines,  to  derive  higher  order 
2 

approximations  of  y^  in  L  (Q) .  However,  it  would  be  difficult  to  justify  a 
better  estimate  than  (3.51).  For  example,  in  order  to  ameliorate  the  counterpart 
of  (3.54),  one  has  to  show  that 


^'•^A  2  -3/2 

1  1  L  (0,T;  H  '  (S))  • 


-»•  0  as  e  ->•  0 


which  is,  by  no  means,  trivial. 

Consequently,  no  attempt  will  be  made  to  pursue  this  any  further. 

3.4.  Hyperbolic  Boundary  Value  Problems 

In  this  section,  an  evolution  problem  of  hyperbolic  type  is  investigated, 
namely,  the  following  boundary  value  problem: 


(y._"  .X)  +  ^  =  (f,x),  (3.76) 
y^(0)  =  h,  h  given  in  V  (3.77) 
y^(0)  =  g,  g  given  in  H  (3.78) 
y^e  L^(0,T;V),  y^S  L^(0,T;H),  f  e  l"(0,T;H)  .  (3.79) 
The  present  analysis  would  be  parallel  to  that  of  Section  3.3.  Hence, 


the  convergence  of  y_  as  c  -»■  0  is  studied.  Then  a  zeroth  order  approximation 


> 


•I  »i 


of  is  constructed,  using  the  weak  limits  of  the  eigenvectors  of  (i.e., 
the  operator  associated  with  the  bilinear  form  +  e  a^(ifi,p) 

(Cf.  Chapter  2).  It  is  well-known  that  hyperbolic  boundary  value  problems 
are  more  complex  than  their  parabolic  counterparts.  Moreover,  not  all 
eigenvectors  of  converge  in  V,  as  proved  in  Section  2.3.  Therefore,  (3.76) 
has  to  be  specialized  to  second  order  operators  A^,  so  that  one  can  specify 
exactly  what  is  needed  for  the  present  analysis  to  hold  true. 

In  the  sequel,  the  following  problem  is  considered 


-2“  ''O  >"£0  "  'o  ™  ‘>0 


—  +  £  Aj  y^l  -  on 


5'cO  ■  °  ^0'  >'£1  ■  "  o” 


^eO  '  ’'el 


y^(0)  =  h(=(hQ,h^)),  h^  €  i=0,l 


-(0)  =  g  on  fi,  g  e  i^(Q) 


e  L'-(0,T;Hj(n)),  €  L"(0,T;L"(rO) 


f  e  L^(0,T;L^(n)) 


where  the  operator  A  is  as  in  Section  3.3. 


Convergence  of  v-  as  e  -*■  0 


This  subsection  commences  with  the  statement  of  the  convergence  theorem 
and  then  a  constructive  proof  of  it  follows  as  in  Section  3.3. 

Theorem  3.5;  Let  y  be  the  solution  of  (3.80)-(3.86) .  Then,  given  a  sequence  of 
e  converging  to  zero. 


Yg.  ^  y  weakly  in  L^(0,T;L^(fi)) 


(3.87) 


’y  ^  9  ? 

^  weakly  in  L  (0,T;L  (Q)) 


(3.88) 


Moreover , 


C  is  a  constant  independent  of  e. 


(3.89) 


k  Ic  k  Ic 

Proof:  Let  }  and  (u  ,'l>  }  be  the  exact  eigenvalue-eigenvector  pairs 

^  k=l  ^  ^  k=l  2 

of  A^,  with  the  eigenvectors  normalized  (in  L  (H))  for  fixed  t.  Using  a  finite 
dimensional  approximation  of  y^  such  as 


k=l  k=l 


(3.90) 


it  is  shown  in  [23]  that 


y™  ->■  y^  strongly  in  L^ (0 ,T;Hq(J^) )  as  m  ->•  +  ® 


(3.91) 


^  .m 

2  2 

— rr  -*■  — y—  strongly  in  L  (0,T;L  (TO)  as  m  -^  + 


(3.92) 


%  _ 


+  v-5' 

^  k=i  ^  ^  k=i  ^  ^ 

{c^}”  ,  {d^}°°  satisfy  the  following  set  of  ordinary  differential  equations 

^  k=l  ^  k=l 

j2  k  . 

k  k  k 

— I  +  yj  =  (f.-i-p 

dt 

c^(0)  =  \  (3 

dck 

-^(0)  =  (8.»^) 


d^d^  ,  . 

- 1  +  d^  =  (f,^^) 

dt2  e  e 

d^(0)  =  (h.v’^) 


d  d  , 

^(0)  =  (g..J)  j 

k=l,2,. . . 

whose  solutions  are  given  respectively  by 

k  k  /  k  ,  k  T  ^ 

c  (t)  =  (h.’j;  )  cos  /u  t  +  — r -  sin  /y  t  +  /  - r -  (f,ip  )dT 

^  ^  0  /U^ 

S  € 


.k.,.  ,,  .k.  /,k.  .  .  ,,k.  .  I  '"c' 

d  (t)  =  (h,*^  )  cos  /X  t  +  -  sin  vX  t  +  - 

e  "  "  Ak  ^  b  Ak 


T  sin  Ak(t-T) 


(f  ,'P^)dT 


Remark  3.22;  Since  f  e  L^(0,T;L^(n)) ,  h  =  (h^.h^)  e  hJcHq)  X  and 

2 

g  €  L  (Q) ,  the  following  inequalities  hold 


iV( 


I 


5:  J  ((f.’O  +  (f.^p  )dt  <  » 

k=l  0 

,  k 

00  <»  ,  A 

S  <  00  ,  I  -r  < 

k=l  "  "  k=l  "  " 


2  ((g.A^  +  (g,^S^)  <  “ 

k=l  ^  ^ 


(3.98) 


(3.99) 


(3.100) 


Now,  let  E  be  a  sequence  converging  to  zero.  It  was  shown  in  Theorem  2.1  that 


■IS  k  k  „,k  ,k  2 

1)  ^  Uq  >  0,  -»■  ijj  weakly  in  L  (Q) 

k  k  k  1 

2)  ->■  0,  linearly  in  e,  -►  </?  strongly  in  11^(0) 


(3.101) 


(3.102) 


from  which  (3.87)  is  deduced. 

Differentiate  (3.93)  with  respect  to  time  and  take  the  limit  as  e  0 

using  (3.101)-(3.102)  to  get  (3.99).  This  limit  is  well-defined  because  of 
(3 . 98)-(3 . 100) ,  which  hold  as  e  ->■  0. 

Consequently,  y  can  be  written  as 


^  k  k  .  _  ,k  ^k 

y=Zc’4;  +  Zdi^’ 


where  {c^}  ,  {d^}  satisfy  (3.96)-(3.97)  after  letting  e  -*■  0,  i.e, 


(3.103) 


V  p  If  r  If 

c  (t)  =  (h,t  )  cos  /p  t  +  -  sin  /p^t  +  /  -  (f  ,'j;^)dT 

0  /  k  0  0  /pi' 

0  0 


d*'(t)  =  (h,¥>*')  +  (g,V5^t  +  /  (t-T)(f,^’')dT 


(3.104) 


(3.105) 


Now  using  (3.93)  and  Theorem  2.1,  one  obtains  (3.89). 
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.2  0 


+  e  A,  y  ,  =  on  0- 
1  si  1  ’J 


(3.107) 


y,i  =  0  on 

0  0  O 

y,i  =  yg  on  R 


y^-lCO)  =  on  Qq 


3^(0)  .  in 


Problem  (3.107)  is  a  nonhomogeneous  boundary  value  problem.  Using  transposition, 

y^,  is  defined  in  a  weak  sense  as  in  Section  3.3. 

•^el 

In  order  to  derive  an  asymptotic  error  estimate  between  y^  and  y^, 
they  need  to  be  redefined  using  transposition.  For  this  purpose,  let 


3X  9^X 

»0  ■  *0  ^  L^q„),  ao  +  .  0  on  I„, 

Jt 


<0-'*  Aq 


^0  ,■> 
—  =  0  on  R,  Xq(x,T)  =  0,  ^^(x,T)  =  0)" 

^0 


(3.108) 


3x  3^X 

<fl  =  (x^:  Xi  e  L^(0,T;H^(fi^)),  €  L^(Q^),  —^+  e  A^X;^  ^  L^(Q^),X^=0  on  Z^UR, 

3t“ 


Xj^(x,T)  =  0,  — ^(x,!)  =  0}^ 


It  can  be  easily  verified  that  y^,  y^  satisfy 


(3.109) 


is  endowed  with  the  topology  carried  over  by  the  mapping 
€  L^(0,T;L^(Qi))  !—  Xi,  i=0,l. 


I  ^  AoXo)dQo  =  /  foXodQQ  -  /  Hq  -^(x,0)df2 


+  /  goXo(x,0)dag  +  e  /  ^Xo  dE.  Vx„e4g 
“n  R  A, 


(3.110) 


f  ^1 


+  /  g^Xi(x,0)dnj^  -  e  /  y^  ^ 


Q  “l-l'-’-'-'-l  -  ^  -cO  ^  dR,  vx^e$ 

g2 

■"  ''oXo>‘'%  '  -  I  \  ^(y-oi^^c 


*  I  SoXo<’‘-°'''“o>  %% 
“o 


(3.111) 


(3.112) 


/■  0  ^1  f  .  ^Xi 

^2  °'-  ^2  **2 


+  /  gj^X2(x,0)dnj^  -  e  /  y 


0  ^^1 

0  3v.  ^X2^'J>2 


R  A, 


(3.113) 


Theorem  3.6;  Let  y^,  y^  be  the  solutions  of  (3. 110)-(3 . 113) .  Then  for 
sufficiently  small  e  the  following  estimate  holds 


ly  -  v^n  -  <  C 
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■  '.s 


a 


Proof :  Subtract  (3.112)  from  (3.110)  and  (3.113)  from  (3.111)  to  get; 


3^ 


I  ~  ^0^^  2  ^  I  3v  ''^O 

Q„  3t  R 


3^X 


0  ’Xl 


I  <^1  -  y“i) (-^  +  £  *iXi)'‘')i  -  -  '  /  (y,o  -  Xo>st 

U-  d  t  R 


dR 


Now  consider  the  following  equations 


ter 

iS- 

> » - . 


3^ 


eO  .  . 

r  "^O^eO  =  ^eO 


3t 


X,o  =  0  on  Eq 


yo  on  Qq 


3X 


eO 


9v , 


=  0  on  R 


Xg.o^x*T)  =  0  on  Qq 


3X 


-(x,T)  =  0  on  Qq 


,2 
c  X 


El  . 

- —  +  e  A  =  y.i 

X^l  =  0  on  E^ 

X^l  =  0  on  R 
X^j^(x,T)  =  0  on 


^1  "^1 


aX 


El. 


"(x,!)  =  0  on 


/ 


(3.115) 


(3.116) 


(3.117) 


(3.118) 


It  can  be  shown  [30]  that 


€  ..  .  i=0,l  . 


(3.119) 


Consequently, 


^eO  „  e  L^(0,T;H^^^(S)) 


(3.120) 


^el  2  -1/2 

—  e  L  (0,T;H 


(3.121) 


From  (3.89),  one  concludes  that 


eO  "2 


<  C, 

1/2  ~  1 


L‘(0,T;H  '  (S)) 


(3.122) 


9y  1 

/T  I  ^  ^  I  2  1/2  —  ^2 

L^(0,T;H 

1  K 


(3.123) 


Let  Xq  =  X^Q  in  (3.115)  and  =  X^j^  in  (3.116)  and  use  (3 . 117)-(3 . 118) , 
(3.122)-(3.123) ,  to  get 


1^0  -  ^  S 

1^.  - 


(3.124) 


(3.125) 


and  hence  one  obtains  (3 . 114)  using  (3 . 124)-(3 . 125) . 


Now  the  weak  limits  of  the  eigenvectors  of  A^,  i.e.,  (3 . 65)-(3 .66)  are 

Ic 

employed  to  solve  (3.112)-(3.113) .  First,  renormalize  {')j  }  .  Then,  the 

k=l 

solution  of  (3. 112) -(3. 113)  may  be  represented  by: 


e  L^(0,T).  Z  /  lc‘'(t)|^dt  < 
k=l  0 


k...  i2. 


(3.126) 


,0  _  _  k  k 

k=l 


d*"  €  L^(0,T),  Z  /  |d‘'(t)|^dt 


k,.,  |2 


k=l  0 
k  k. 

In  order  to  obtain  c  (t) ,  d  (t),  let 

e 


(3.127) 


XgCx.t)  =  9(t);(,Q(x),  eec^([o,T]),  0(1)  =  ^(t)  =  0 
Xj^(x,t)  =  v(t)v’^(x),  v€C^([0,T]),  v(T)  =  =  0 


in  (3.112)-(3.113)  to  get 
T  2  T 

/  =  /  (fQ,'i;J)0dt  -  (hQ,'^Q)  1^(0)  +  (gQ,'^J)0(O) 


(3.128) 


0  Aj£v)dt  =  /  (f^,^J)9dt  -  (h^.^J) 


+  (8,  ,>^1  )'j(0)  -  e  /  (  Z  c  /  -r— ^  dS)  vdt 

0  ^=l  S 


(3.129) 


which  are  equivalent  to 


d^c*"  ^  k  k  k. 

^  +  “o  -=  ■  «0’'V 


c’'(0)  -  (h.,*!;) 


(3.130) 


dt 


dNo)  =  (h.,>^h 

e  1  i 


d  d 

^(0)  =  (g^..p 


(3.131) 


Remark  3.24;  Remarks  (3. 16)-(3.17)  also  apply  in  this  case. 


Remark  3.25:  If  h  =  (h^.h^)  ^  belongs  to  some  technical 

difficulties  would  be  encountered  in  defining  (3.106). 

Remark  3.26:  The  difference  between  the  estimates  (3.63)  and  (3.129)  is 
due  to  the  fact  that  the  solutions  of  (3.63)-(3.64)  are  more  regular  than 
the  solutions  of  (3 . 117)-(3 . 118) .  For  further  inquiry,  the  interested  reader 
is  referred  to  [23,30].  ° 


An  example  is  now  given  to  illustrate  the  computational  aspects  of  the 
approximation  of  y^.  This  example  is  also  considered  in  applications  in 
Chapter  5. 


Example  3.7:  Let  f„=l,  f,  =0,h=0,g=0.  Then  (3 . 80) -(3 . 84)  become 
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3t 


^'>'1 

ot 


3x‘ 


f-  =  1  on  (-l,0)X(O,T) 


3^ 


-  e 


cl 


=  0  on  (0,1)X(0,T) 


3x‘ 


y^Q(0,t)  =  y^^(0,t) 


"^cl 

-^(0,t)  =  c  -^(O.t) 

OX 

y.Q(x,0)  =  y^^(x,0)  =  0 

9y„n  ^y-1 

-^(x,0)  =  -^(X,0)  =  0 


Recall  that 


=  ((2  k-1) 


1)^  =  'll  cos  /UqX 


=  (kTi) 


^ ^  =  /l  sin 


The  solutions  to  (3 . 130) -(3 . 131)  are  readily  calculated  and  are  given  by: 
=  "V  IT  3/2 

(’.q) 


d^(t)  =  ^  (1  -  cos  X^t)  S 

’'•*1 


+  2/2  A^c  T 


?.-l 


l\l2  1,  k  .  ''^1'-’^  • 


^=1  (m3)-"“(Uo-v,c) 


(3.132) 


Therefore,  y  is  written  as 
e 


.k-1 


°  ^  ^  (1  -  cos  /yjt)  cos  /u^x 


^0  ;  k,3/2 

k=l  (Uq) 


(3.133) 


0  ,  .  k  {  ,.  _ 


4,-1 


^  k=i  V  -  !i=i 


(3.134) 


+  e  Z 


(-1) 


£,3/2,  £  ,k  , 


-  cos 


£=1  (Uo)^'^(vX,c) 


/A^t) 


sin  A^x 


Remark  3.27;  The  counterparts  of  Remarks  3.19-3.20  are  also  applicable  in  this 
example . 


3.5.  Concluding  Remarks 

In  this  chapter,  the  convergence  of  the  solution  of  three  classical 

boundary  value  problems  (namely  elliptic,  parabolic  and  hyperbolic)  as  e  -»■  0 

has  been  analyzed,  using  the  spectral  analysis  undertaken  in  Chapter  2.  For 

elliptic  boundary  value  problems,  it  is  found,  that  by  modifying  the  weak  limits 

appropriately,  a  strong  Laurent  series  expansion  of  y^  can  be  derived.  For 

parabolic  and  hyperbolic  boundary  value  problems,  zeroth  order  approximations 
2  2 

in  L  (0,T;L  (fi) )  were  easily  constructed,  using  only  the  weak  limits  of  the 
eigenvectors  of  A^. 

Several  examples  were  solved  to  illustrate  and  clarify  the  aspects 
of  the  problems  at  hand.  Table  3.1  summarizes  the  properties  of  boundary  value 


problems  investigated  in  this  chapter,  in  their  simplest  form. 


SUMMARY  OF  BOUNDARY  VALUE  PROBLEMS  INVESTIGATED  IN  THIS  CHAPTER. 


Zeroth  order  approximation 
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The  approach  followed  in  the  present  exposition  is  general  enough 
to  encompass  many  stiff  operators.  As  an  example,  many  boundary  value  problems 
involving  the  stiff  operators  considered  in  Sections  2. 3-2. 5  can  be  approximated 
using  the  same  concepts  developed  herein. 


CHAPTER  4 


SUBOPTIMAL  CONTROL  OF  STIFF  SYSTEMS 


4.1.  Introduction 

In  Chapter  2  the  spectral  decomposition  of  a  class  of  stiff 
operators  A^ ,  including  the  convergence  of  their  eigenvalue-eigenvector 
pairs  as  e^O,  has  been  analyzed.  Using  the  results  of  this  investigation, 
some  classical  boundary  value  problems  involving  the  aforementioned 
operators  were  studied  in  Chapter  3.  The  convergence  of  their  solutions  as 
e 0  were  analyzed.  Then,  asymptotic  approximations  of  these  solutions  were 
constructed,  using  the  weak  limits  of  the  eigenvectors  of  A  as  e-^0. 

'C 

Asymptotic  error  estimates  were  also  obtained. 

In  the  present  chapter,  some  control  problems  with  quadratic  cost 
functionals  are  considered.  The  results  derived  in  the  previous  two 
chapters  are  used  to  investigate  these  problems.  The  objectives  of  this 
chapter  are: 

1.  to  obtain  information  about  the  behavior  of  the  optimality  system  as 

2.  to  "approximate"  the  state  and  the  control  of  the  system  for  small  values 
of  e . 

The  control  problem  of  distributed  parameter  systems  is  formulated 
in  many  books  and  manuscripts  such  as  [2,8,23-30],  to  name  a  few.  In  [8], 
a  semi-group  approach  is  followed.  However,  a  variational  approach  is 
chosen  in  [23]  and  the  subsequent  references.  Recently,  several  results 
about  Dirichlet  boundary  control  in  parabolic  and  hyperbolic  systems  have 
appeared  in  the  literature  [7,20,21].  However,  the  assumptions  made 
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therein, concerning  the  coefficients  of  the  operator  A^,  are  more  restrictive. 
Hence,  a  control  problem  with  Neumann  boundary  control  is  considered  in  the 
sequel. 

The  formulation  of  [23]  seems  to  be  adequate  for  the  presentation 
herein.  Consequently,  the  control  problems  to  be  considered  are  adapted  from 
there.  In  this  chapter,  the  control  of  a  class  of  parabolic  systems 
is  investigated.  Two  types  of  control  are  considered.  In  Section  4.2, 
the  control  is  distributed.  In  Section  4.3,  the  control  is  of  Neumann 
type,  exercised  through  the  boundary.  In  each  section,  the  problem  formu¬ 
lation  is  first  presented.  Then  the  convergence  of  the  state  and  the  costate 
as  e -► 0  is  studied.  Their  asymptotic  approximations  are  then  constructed, 
using  the  approach  developed  in  Section  3.3.  In  Section  4.4,  some  concluding 
remarks  are  given. 


4.2.  A  Parabolic  Problem  with  Distributed  Control 
4.2.1.  Problem  formulation 

Let  H  =  L^(^^),  V  =  (fi)  where  U  SC  ,  with  boundary 

r  *  FgU  as  depicted  in  Figure  2.1.  Let 


14  44 

aui'P,'.!)  -  I  1  I 
^  i=l  j=l  ^ 


<*> 


(4.1) 


where  a^j (x)  satisfy  the  conditions  of  Remark  2.9. 
control  problem 


Consider  the  following 
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inf  J  (V)  =  j  lly  (T)-zJ|^  +  j  lly  -z 
veL2(0.T;H)  ^  ^  2  II  e  dll^2 

+  j  (Nv,v)  2 

L  (O.T;H) 


(4.2) 


subject  to 

•  (y^.s5)  +  agCy^.ss)  +  ea^(y^,.^)  =  V^ev  (4.3) 

•  y^(0)  =•  h,  h  and  are  given  H  (4.4) 

2 

•  f  and  are  given  in  L  (0,T;H)  (4.5) 

•  y^eL^(0,T;V)  (4.6) 

2  2 

•  N  is  a  given  operator  in  JC  (L  (0,T;H);  L  (0,T;H)),  which 

is  herraitian  and  positive  definite.  (4.7) 


Under  these  assumptions,  the  above  control  problem  admits  a 

2  2 

unique  optimal  solution  [23]  {y^,u^}€L  (0,T;V)^  L  (0,T;H)  for  fixed  e, 
characterized  by  the  following  optimality  system 

(y^.^)  +  aQ(y^,sj)  +  eaj^(y^,'^)  =  (f,i^)  +  (u^,'^), 

(-p^,v)  +  a^(p^,-v)  +  ca^(p^,l)  =  (y^-z^,..), 
y^(0)  -  h,  pjT)  =  y^(T)  -  z^ 
u  *  -N~^p 
y^.P^e  L^(0,T;V) . 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


The  optimality  system  (4. 7) -(4. 10)  can  be  decoupled  through  the  affine  map 


where  P  ,  r  are  defined  by 

t  e 


-P'+PA  +AP  +P  N'^P  =  I 

e  e  c  e  e  e  e 


P  (T)  -  I 

e 


-r'  +  A  r  +  P  n“  r  =  P  f  - 

e  c  e  e  e  e  d 


where  A  is  the  operator  associated  with  the  bilinear  form  a  = 

e  G 

,'li) .  Consequently,  the  optimal  control  given  by  (4.10) 
can  be  written  in  the  feedback  form  as 


u  *  -N  (P  y  +r  ). 
e  e  e  e 


The  properties  of  P^  and  r^  are  summarized  in  [23],  Theorem  4.4 
p.  148,  some  of  which  are 

P^(t)€i:(H;H)  (4 


P*(t)  =  P  (t),  (P  (t)s,s)  >  0,  VsEH 

e  e  G 


r^eL  (0,T;V).  (4 

r  Ic  Ic  ^ 

Now  let  »Xj.  the  eigenvalue-eigenvector  pairs  of  the  operator  A_ 


we  V 


such  that 


(x^fX'^)  =  '5^'^  (Kronecker  delta), 
e  c 


(4.20 


For  simplicity,  assume  N  =  ol,  p  >  0  and  the  eigenvalues  of  are  not 


repeated.  Then  the  operator  P  (t)  can  be  expressed  as 

e 


P  (t)VJ  =  E  E  P^^X^C'^.xJ'). 

e  i=l  j=l  e  e  ^ 


V'i’G  H 


where  {P  .  .  ,  satisfy 
e  i.J=l 


dp 


ij 


e  .  /  i  .  ix  ij  .  1  ..  ki  kj 
■T—  +  (y  +Y'^)p  +  —  E  p  p 

dt  etc  p  1^=1  £  e 


=  5^^ 


pJ^(T)  =  6^^ 


i.j  =  1.2,. 


It  can  easily  be  verified  [8]  in  this  case,  that 


p^^  =  0  if  i#j 

c 


is  a  solution  of  (4.22).  Therefore, (4 . 22)  reduces  to 


ii 


^  1  11  .  1  ,  11^2 
-7- —  +  2y  p  +  —  (p  ) 
dt  p  G 


=  1 


p“(i) 


i  *  1,2, .. . 


The  function  r^  can  be  expressed  as 
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r  =  -  s  X 

c  i=i  G  G 


(4.21 


(4.22 


(4.23 


(4.24 


(4.25 


where  {s  }.  ,  satisfy 
£1=1 


“"e  .  .  ii,  i 

-rr-  +  (y  +—  P  )s 

at  epee 


i.  ii  .  i, 
(f.xpp^  -(Zd.X^) 


s^T)  =  -(Zf.xJ) 


i  =  1.2,, 


4.2.2.  Convereence  of  the  state  and  the  costate  as  e ->■  0 


The  convergence  of  and  as  e  ^  0  is  summarized  in 
Theorem  4.1;  Let  y^  and  be  the  solution  of  the  optimality  system  (4.8) 
(4.11).  Then  as  e->0, 

2 

y^  ->■  y  weakly  in  L  (0,T;H) 


p^  ^  p  weakly  in  L  (0,T;H) 


L  (0,T;V) 


L  (0,T;V) 


J^(u^)  ^  J(u) 


where  y  and  p  satisfy 


(y' ,'fi)  +  Si^iy,^)  +-  (p,v?)  =  (f.v?). 


(-P'.  l')  +  agCp,’!')  -  (y.i')  =  -(Zd»  J-)* 


v^ev 


V 


y(0)  =  h. 


p(T)  =  y(T)  -  Z-. 
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Proof:  Let 


”  k  k  ^  V  ^k  k 
y  =  Ec';  +  I  d 

^  k=l  ^  e  k=l  ^  ^ 


;  k  k  “  .  k  k 
p  —  Za'4/  + 

e  k=l  e  e  k=l  e  e 


where  „d  satisfy 


dc 


e.kk^lk  /fk, 
dT  •"  “s't  *S\’ 


da 


e  ,  k  k  k  /  k, 
■j—  +  M  a  -c  =  -(Zj.ijy  ) 
dt  £  e  e  d 


c^(0)  »  (h,;i,J),  a^(T)  =  c^(T)  - 


1 


_ e  .  ,k,k  ,  1  ,  k  ek, 

-r—  +  Xd  +;rb  =  (f  ,'f  ) 
dt  G  e  N  e  g 


dt  £  G  G  d  £ 


dNo)  =  (h.^S,  bNx)  =  d^(T)-(z.,^S  1 

G  G  G  G  t  G  ' 


k  =  1,2,... 


The  above  coupled  equations  have  unique  solutions  for  any  value  of  g.  As 
£-*■0,  using  the  results  of  Theorem  2.1,  take  the  limit  in  the  above  equations 
to  derive  easily  the  conclusions  of  Theorem  4.1. 

For  a  more  general  approach  that  is  applicable  when  N  is  a 
function  of  x  or  the  control  is  exercised  on  the  boundary,  consult  [24]. 

As  a  result  of  this,  the  Riccati  operator  P^(t)  and  the  function  r  (t)  also 
converge  in  some  sense. 


Theorem  4.2:  Let  P  (t) ,  r  (t)  be  the  solutions  of  (4. 13)-(4. 14) .  Then  as 
-  e  e 

e  -*•  0, 


P  (t)  P(t)  in  £,  (H;H) 

e 


r^(t)  -►  r(t)  weakly  in  L  (0,T;H) 


P(t),  r(t)  can  be  written  as 


P(t)  = 


'Po(t)^ 

,Pl(t)y 


r(t)  = 


\(t)> 

.ri(t)y 


where 


Pf,(t)x  =  2  p^^v^Cx./) . 

u  X=1  ^ 


vxeH 


T,  /  x  _  li  i,  ix 

P,  (t)x  =  P,  (.Xi'f  )  . 


VxeH 


r  (t)  =  I 

1  i=l  1 


^^k^i=r  ^®k^i=r  satisfy 


A 

^0  ^  ,  i  ii  ^  1  ,  iix2 

-dT-  ^  +  -  (Pq  ) 


p:^(t)  =  1 
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=  0  on  I 
0  0 


3y 


5^ 


•  =  0, 


Sv. 


=  0  on  R 


■  %■  '  "fO  “o 


^  0 


^  0 

.  .  0 


.  1  0 
^1  "n  Pel 


^el-'d 


on  Q, 


on  Q, 


0,  P*^  = 

0 

on 

“1 

0 

0  0 

0 

* 

^0’  Pel  “ 

Pq 

on 

R 

/ 

y,Vo) 

\ 

on 

"l 

\ 

P,^(T) 

^fl 

on 

"l 

I 


(4.41, 


(4.42] 


(4.43) 


(4.44'' 


(4.45 


where  the  solution  of  (4.43)-(4.45)  is  defined  as  in  Section  3.3.2  by 
transposition. 

0  0 

Theorem  4.3;  Let  (yg.»Pg.^  the  srlution  of  (4.8)-(4.10)  and  (y^»P  1  be 
the  solution  of  (4.40)-(4.45) .  Then  the  following  estimates  hold  for 


0  <  e  <<  1 


el 


(t)  =  E 


i=l  el 


(4.f 


with  respectively, 


i  i” 


dp 


11 


el  .  o,i  ii  ,  1  /  ii.2  , 
dT  +  +  p  -* 


(4.: 


p“(i)  -  1 


ds  i 

el  ,  .1  il-N  i  .is  r  0  3<fi  ii  \ 

■"  Ppp^pi '  «*>'  ’  I  ^0  SET  " 

S  A- 


dt 


el 


-e  /  p’ 


0  3^ 


S  0 


dS 


sJ^(T)  =  -(z^.^S. 


(4.5 


The  presence  of  and  p^  in  (4.56)  implies  that  the  optimality  systems 


given  by  (4.40)-(4.42)  and  (4.43)-(4.45)  have  to  be  solved  sequentially. 


0 


However,  it  is  clear  from  the  analysis  of  Section  3.3  that  setting  y^  and 


p2  to  zero  in  (4.44)  induces  errors  in  y^,  and  p^ ,  no  more  than  O(’^'e)  . 
0  •^el  *^cl 


Consequently,  it  may  be  desirable  to  set  them  to  zero.  Another  possibility 


may  be  that  no  control  is  exercised  on  In  such  a  case,  it  may  be 


/  ^  CO 

advantageous  to  synthesize  a  feedback  law  of  the  form  (4.52),  where 


0 


are  computed  from  (4.56)  with  p^  set  to  zero. 


Example  4.1:  (Cf.  Example  3.6)  In  this  example,  the  following  control 
problem  is  analyzed: 


ini 


v£L^Q)  «)  I-  «) 


(4.5 
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subject  to 


’>'£0 


0  t 


cfX 


o  =  V 
2  0 


3y 


el 


3^y 


el 


3t 


9x 


on  Q 


0 


=  on  Q 


1 


(4.58) 


=  yg-i(o.t) 


=  0 


(4.59) 


3x 


3y  1 

(O.t)  =  e  (0,t) 


(4.60) 


=  U  y^^(x,0)  =  1. 


(4.61) 


A  suboptiraal  feedback  control  law  as  outlined  previously  would  be  given  by 


0  1  „  0 
"0  ■  -  F 


^  1  =  —  p  ly  1  • 

el  0  el  el 


) 


(4.62) 


Recall  that 


=  ((2k-l)-2)  ,  0  “  cos(2k-l)^  x 


\\  =  (kTT)^. 


'f  ^  =  /2  sinkiTx. 


Hence,  (A. 62)  can  be  rewritten  as 


0  _ 
“o  ~ 


el 


1  ^  i,  0  i.  i 

S' Ji^oV-V- 0 


1  "  i  ,  9  i.  i 


p  i=i  - 


,  OD 

Where  ^Pel^i=l 


^^0  .  „  i  i  ^  1  i  ? 
■dT  ■"  ^^0^0  ■"  p 


=  1 


Pq(T)  =  0 


dp^ 

£l  .  .1  ,  i  .2 

*7  '‘’eI* 

p‘l<T)  -  0. 


=  1 


In  this  case,  snd  be  computed  in  closed  form 


1 


1  i  -°o^b-T) 

1  Vo<^  -‘> 


Pq-  » 


^0  ®o 


i  .i 


Pel  “  P 


B  ,A  , (e 
el  el 


-1) 


,1  -\i'^-''>  .1 

-hi 


as 


where 


130 


.i  _  i  /  iv2  ,  1 

0  0  0  P 


_  i  /.  i.  2  1 

Bq-  -Uo+  /(Pq)  +- 


(4.68) 


el 


=  -xj-e-  /(\^£)^+^ 

i  1  P 


=  -X 

cl 


te  +  /(X^'^  •  ^ 


(4.69) 


The  numerical  results  are  discussed  in  Chapter  5  for  various  values  of  c 
and  T. 


4 


r 


.3.  A  Parabolic  Problem  with  Neumann  Boundary  Control  and  Boundary 
Observation 

4.3.1.  Problem  formulation 

Let  H  =  L^(rO»  V  =  H^(n)  where  SCir’^  with  smooth  boundary 

and  interface  S  as  indicated  in  Figure  2.1b.  Let 


where 


A  = 

(T 


0 

eA 


n  n 
A,  =  Z  E 
^  i=l  j=I  iJ 


3  k  .  .  3 

a. , (x) 


3x.  ’ 

J 


k  =  0,1 


Let 


k 

"ij 


satisfy  the  conditions  of  Remark  2.9. 


=  n^x(o,T),  i  =  0,l 

Z  =  rx(o,T) 


R  =  Sx(0,T) . 


Now  consider  the  following  control  problem 


inf  J  (v)  =  / 

veL2(Z)  ^ 


dZ  +  (Nv,v) 


subject  to 

•  (y'.'^)  +  a  (y  ,'P)  +  ea  (y  ,</»)  =  (f  .v^)  +  (v.V’)  ,  W€  V 

e  u  e  1  e 

•  y  (0)  *  h,  h  given  in  L^(Q) 

•  e 

•  f  =  (f^,fp  .  given  in  L^(0,T;H)  such  that 

(f.,1)  ,  =  0,  i  =  0,l 

L  (n^) 

•  given  in  L^(Z) 

2  2 

•  N  is  a  given  operator  in  £ (L  (Z) ;  L  (Z)), 
which  is  hermitian  and  positive  definite. 


(4.70) 


(4.71) 

(4.72) 


(4.73) 


(4.74) 


(4.75) 


Under  these  assumptions,  the  above  control  problem  admits  a  unique  optimal 

2  2 

solution  {y^,u_}€L  (0,T;H)  ’<  L  (Z)  for  fixed  e  [23],  characterized  by  the 


following  optimality  system 


(y  +  3l  Ay  +  ea  (y  ,^)  =  (f,<;^)  +  (u  ,^)  , 

^  UE  LE  E_<-/„\ 

L  (i  ) 


Vw-  £  V 


>  (A. 7 


"(p'.A)  +  a  (p  ,v)  +  ea/p  ,v)  =  (y  -z  ,v) 
e  0  c  1  E  Ed 


Vv-  V 


y^(0)  =  h,  p^(T)  =  0 


y^.p^GL  (0,T;V) 


(A. 7 
(A. 7 


where 


u  =  -N  ^p  j 
e  e 


(A.  7 


This  optimality  system  can  be  decoupled  as  in  Section  A. 2  by  the  affine  map 


p  =Py  +r,  P*  =  P 

e  e  e  e 


(A. 81 


where  P^  and  r^  satisfy,  respectively. 


(-p'^'P  ,v)  +  »P^r)  +  ea^Cv’  +aQ(P^i^>  ,v)  +  ea^(P^i^  ,y) 


+  (n'^p  ,p  ■;.)  ,  =  ,  ¥^,vev 

^  L^(r)  L^(r) 


>  (A. 8 


P  (T)  =  0 
e 


(-r'  +a-(r  ,v5) +ca,  (r  .VJ)  =  (P  f,i^)-(z,  ,v5)  „ 

E  0  £  1  £  E  d 


Vv?€  V 


)  (A. 8: 


r  (T)  =  0. 

e 


Ic  k 

Now  let  be  the  eigenvalue-eigenvector  pairs  of  ,  i.e., 


(X^.v')  +  ea,  =Y^(x^.v"),  (x^tX^)  =  <5^^  , 


Vv'S  V, 


(A. 8 


For  simplicity,  let  N*pl,  p>0,  then  P  and  r  can  be  written  as 

e  e 


P  f  =  Z  E  .X^)  . 

^  i=l  j=l  e  £  e 


Vy€  V 


r  =  E  s\^ 
e  i=l  e  e 


ij 


where  .  ,  and  {s^}.  ,  satisfy 

■^e  i,j=l  e  1=1 


4.  /  i.  V  ^  ki  ilj  r  k  Jl  r  i  j  ^ 

+  +Y Jp^  I  P^  /  X  X  dr  =  /  X 


dt 


e  ^  P  k=l  £=1  £  £  f  £  £ 


e  £ 


Mr 


> 


pJ^(T)  =  0. 


l,j  =  1,2,. 


where 


Pe  ’ 


Vi.j 


ds  .  , 

e  ,  1  i 

~JT~  +  Y  s 
dt  e  e 


'T  ^ ^  /  i\ 

^  P-  (f  .X_)  -  o 

i=l  e  e  d  e 


(4 


s^^T)  =  0. 

e 


Remark  4.2;  Since  V  =  (fi)  ,  the  operator  has  an  eigenvalue  which  is 
zero  ind  its  corresponding  eigenvector  is  a  constant  on  n.  By  assumption 
(Cf.,  (4.73)),  they  are  excluded  from  (4.83). 

Remark  4.3;  Observe  that  (4.86)  is  not  decoupled  as  in  Section  4.2. 


It  is  straightforward  to  prove,  as  in  Section  4.2,  the  following 

theorem. 

Theorem  4.4;  Let  and  be  the  solution  of  the  optimality  system 
(4. 76) '(4. 79)  .  Then  as 

2 

y^  y  weakly  in  L  (0,T;H) 

2 

p^  -»■  p  weakly  in  L  (0,T;H) 

J  (u  )  ^  J(u) 

e  e 

where  y,p  satisfy 

(y',',i>)  +  a^(y,<^)  +  ^  (p.i^)  „  =  (f,'^), 

lUt) 

-(p'.v)  +  a  (p.'iO  -  (y,v)  «  =  „  ,  V|€V 

L^(r)  L^(r) 

y(0)  =  h.  p(T)  =  0. 

Proof :  See  proof  of  Theorem  4.1. 

Remark  4.4;  As  in  Section  4.2, 

P^(t)  *  P(t)  in  I  (H;H) 

2 

r_(t)  --  r(t)  weakly  in  L  (0,T;H) 

where  the  pair  <P,rr  satisfy  (4.8l)-(4.82) ,  respectively,  after  letting 
£-►0.  Note  that  P.r}  decouples  the  optimality  system  limit,  i.e.,  (3.88). 
Furthermore,  P(t)  and  r(t)  may  split,  as  in  Section  4 . 2,  in  some  instances 
(depending  on  the  forcing  term  f) ,  because  the  influence  of  the  null  eigen 
value  of  is  excluded  by  assumption  (Cf.,  Remark  4.2). 


Asymptotic  approximation  of  y  and 


Let  “  ^^O’^eP  ‘lenote  the  zeroth  approxima 

2 

tions  in  L  (0,T;H)  of  y^  and  p^,  respectively.  For  0 < e <<  1,  they  are 
defined  by 


0  4.  A  0  c 

TT  Vo  =  ^0  Q( 


0 

■  IF  Vo  “  °  on  a 


21^  0 
3v  '^O’ 

"^0 


^o-'d 


on  E 


yo(0)  =  hg. 


Pq(0)  =0  on 


^0  =  -  N  Po 


+  A  0  f 

Tr+  Q] 


—  °  '5] 


0  0 


0  0 


=  ^0’  Pel  “  Pq  ^ 


y°^(0)  =  h^.  p°^(T)  =  0  on  r,^ 


where  the  solution  of  (4.93)-(4.95)  is  defined  using  transposition. 


Theorem  4.5;  Let  be  the  solution  of  (4 . 76) -(4 . 79)  and  ^yg.»Pg.^ 

be  the  solution  of  (4.89)-(4.95) .  Then  the  following  estimates  hold  for 
0  <  e  «  1 . 


2  - 
C.  b  -  ^  m  A 


L  (0,T;H) 


L  (O.T;H) 


lJ^(u^)-J(uo)|  <  036 


Cj^, 03,03  are  some  constants  independent  of  e. 

Proof ;  See  proof  of  Theorem  4.3. 

As  before,  (4 .89) -(4 .95)  can  be  solved,  using  the  weak,  limits  of 


the  eigenvectors  of  the  operator  .  It  can  be  easily  shown  as  in  Section 

Ic  Ic 

2.5.1  that  the  weak  limits  of  decomposed  into 

,  whose  weak  limits  satisfy 
e  e  k»0  e  e  k»l  ^ 


V’  Q  =  constant  on 


A*^^  -  X*'/ 

V  1  11 


k  =  0,1,2, 


(4.100) 


k  =  1,2,...  ^ 

where  the  null  eigenvalue  and  its  corresponding  eigenvector  are  included 

k  1 

in  (4.99).  The  constant  in  (4.99)  is  chosen  such  that  'f  €h  (H)  ,  k=0,l,2. 
However,  because  of  (4.73),  =  0  and,  therefore,  (4.99)  has  to  be 

modified  into 


”0  ■  “ 


.  -k  ,k.k 

Vi  ■  h'l 


on  n. 


on  Q, 


(4.101) 


k  =  1,2,.. 


Now  (4. 89) -(4. 91)  can  be  decoupled  by  the  following  map 

0  D  0  ... 

Po  ■  *  ''0 


with  Pq  and  r^  written  as 


00  00  i*  ^  . 

0  1=1  j  =  l 


s. 


(4.102) 


(4.103) 


(4.104) 


where  .  ,  and  {s„}.  ,  satisfy 

*^0  i,j  =  l  0  1=1  ^ 


dp 


CO  oo 


_  ki  c  k 

dt  ■  '"o'  ■"  p  k=l  l^l^O  '"O  '■  ^ 


r 


/  v\^dr 

r 


Pq^(T)  =  0. 


Pq  =Po 


i.j  = 1,2, ., . 


(4.105) 


ds 


dt 


0  i  i  ”  ii/^  Us  /  is 


L  (r) 


(4.106) 


=  0«>  -  0 


i  =  1,2,... 


Remark  4.5;  Note  that  (4.93)-(4.95)  are  decoupled.  Hence  there  is  no 
0 


need  to  compute  p^  for  control  purposes.  As  before,  (4.89)-(4.92)  and 

(4.93)-(4.95)  have  to  be  calculated  sequentially  in  time  because  of  (4.94) 

0 

However,  by  setting  =  0  ori  R,  a  cost  no  larger  than  0(v'£)  is  incurred 
since  the  solution  of  (4.93)-(4.95)  does  not  influence  (4.92).  Therefore, 
it  is  rational  to  do  exactly  that  to  avoid  the  sequential  computation. 

=  0,  then,  by 


Remark  4.6:  If  the  condition  (4./3)  were  read  (f,l)  , 

1  ,  L^(Q) 


letting  f 


0  meas .  il 


0 


meas.  .1, 


in  (4.71),  some  "controllability 


difficulty  would  be  encountered. 


Example  4.2:  (Cf.,  Example  3.6)  Let 


sin  2:tx  for  x€(-l,0) 


Zd  =  0,  f  =  < 


h  =  0 


for  X  £  (0 , 1) . 


Then  the  control  problem  (4. 70) -(4 . 77)  becomes 


int  J_(v)  =  /  (y  (-l,t))  dt  +  p  /  (v(t))  dt 
v€L  (0,T)  ^  ^ 


subject  to 


3y  „  3  y  rv 

— ^ - “  ^0  (-l,0)x(0,T) 


3t 


3x 


=  °  (0,l)x(0,T)  j 

3x 


y^QCO.t)  =  y^^(0,t) 


3y  3y  , 

(O.t)  =  e  (0,t)  ) 


3X 


3y 


rO 


dX 


(-l,t)  =  V 


y^(x,0)  =  0. 


A  suboptimal  feedback  control  law  can  be  synthesized  as 


0  1  /T,  0  .  ^  I 

“O  "  ■  p  ^^0^0 I 


x=-l 


where  Pq  and  Tq  are  given  by  (4. 103)-(4. 104)  where  {v  are 

r“  ri-i® 

{v2  cos  i-Tx  } .  ,  and  IPn  r.  .  i»  ,  satisfy 

1=1  i,j  =  l  0  1=1 


^^0  .  ^2  ^  ,2,  ij  ^  2  !!  ”  ,  ,,k+i  ki  ;j 

-T—  +  ((it)  +  (j  ^)  )p^  T  ^  P  P 

dt  U  p  1^=1  .1=1  0  0 


=  (-1) 


i+j 


(4.107) 

(4.108) 

(4.109) 

(4.110) 

(4.111) 

(4.112) 


(4.113) 
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ds;  ,  . 

0  ,  , .  .2  1 

^  ‘o 


s  (T)  =  0. 
0 


!!  ±%  2(l-(-l)  ) 

^=1  TT(3i^-4) 


(4.114) 


Remark  4.7:  The  present  methodology  for  approximating  control  problems 
cannot  be  applied  blindly  because  it  may  yield  erroneous  results.  Many 
factors  can  limit  its  applicability,  depending  on  the  control  problem  at 
hand.  Two  prominent  factors,  which  are  central  in  any  control  problem,  are 
the  type  of  control  and  the  type  of  observation.  One  problem  where  this 
methodology  fails  is  the  elliptic  stiff  control  problem  with  Neumann 
boundary  control  and  Dirichlet  boundary  observation  considered  in  [29], 
page  323. 


4.4.  Conclusion 

In  this  chapter,  the  control  of  two  stiff  systems  was  considered. 
Using  the  concepts  developed  in  the  previous  chapter,  suboptimal  feedback 
control  laws  were  derived  for  these  problems  for  small  values  of  the  para¬ 
meter  £. 

It  was  shown,  that  the  approximations  of  the  state  and  the  costate 
are  easy  to  obtain,  provided  some  care  is  taken,  depending  upon  the  specific 
problem  at  hand.  It  is  safe  to  claim,  based  on  the  present  results,  that 
these  approximations  alleviate  stiffness  for  most  problems  with  meaningful 
disturbances.  Control  problems  with  Dirichlet  boundary  control  are  more 
complex.  For  practical  classes  of  control  inputs,  the  state  is  not 


"sufficiently"  regular  [22,23].  Therefore,  the  control  space  has  to  be 
restricted  to  obtain  meaningful  results.  As  a  consequence,  the  feedback 
synthesis  of  the  control  is  involved.  One  remedy  for  this  dilemma  is  to 
assume  that  the  coefficients  of  the  operator  A^,  as  well  as  the  boundary 
where  the  control  is  exercised,  are  more  regular.  In  turn,  this  assumption 
restricts  the  number  of  problems  that  can  be  considered.  For  example,  a 
possible  class  of  stiff  control  problems,  with  Dirichlet  boundary  control 
that  may  be  investigated, is  the  class  of  problems  when  the  operator  A^  is 
as  given  in  Section  3.5.3. 


CHAPTER  5 


APPLICATIONS  TO  HEAT  TRANSFER  AND  ELECTROMAGNETICS 

5.1.  Introduction 

There  are  numerous  dynamical  systems  whose  evolution  can  be 
modeled  better  by  partial  differential  operators  such  as  heat  transfer, 
electromagnetic  wave  propagation,  chemical  processes,  elasticity,  just  to 
name  a  few. 

The  introduction  of  a  (or  several)  small  parameter  e  may  have  a 
physical  meaning  of  a  small  conduction  (or  convection)  coefficient  in  heat 
transfer  or  a  small  permittivity  in  electromagnetics.  It  may  also  be 

completely  artificial,  such  as  in  penalized  and  regularized  problems. 

In  the  previous  chapters,  the  theoretical  implications  of  letting 
0  in  some  of  these  models  have  been  studied.  In  the  present  chapter,  two 
specific  examples  of  such  models  are  considered.  The  first  example  describes 
the  heat  conduction  in  a  one-dimensional  rod,  made  of  two  interfaced  media, 
having  heat  conduction  coefficients  of  0(1)  and  0(e),  respectively.  The 
second  example  considers  the  propagation  of  an  electric  field  in  a  one- 
dimensional  waveguide,  consisting  of  two  interfaced  media,  having  permit¬ 
tivities  of  0(1)  and  0(e),  respectively.  Most  of  the  interpretations  given 
in  the  sequel  are  of  general  nature  and  hence  applicable  in  many  other 
related  problems. 

This  chapter  is  organized  as  follows.  In  Section  5.2,  physical  inter 
pretations  of  the  results  obtained  in  Sections  2.2  and  2.5  concerning  the 
convergence  of  the  eigenvalue-eigenvector  pairs  of  stiff  operators  are  given. 
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In  Section  5.3,  the  asymptotic  approximations  of  the  solution  of  the  boundary 
value  problems  of  Examples  3. 6-3. 7  are  compared  with  finite-dimensional 
approximations  of  these  problems  for  different  values  of  s.  In  Section  5.4, 
the  control  problem  of  Example  4.1  is  solved  numerically.  The  last  section 
contains  some  concluding  remarks. 


5.2.  Physical  Interpretation  of  the  Limits  of  the  Eigenvectors  of 
Stiff  Operators 

Physical  interpretations  of  the  convergence  of  the  eigenvectors 
of  stiff  operators  as  e^O  are  given  within  the  framework  of  the  examples 
discussed  in  Chapter  2. 

The  operator  in  Example  2.1  may  represent  the  heat  diffusion  • 
in  a  slab  occupying  the  space  HCr”  (n  <  3)  ,  composed  of  two  interfaced  media 
having  dif fusivities  0(1)  and  0(c),  respectively.  The  complement  of  the  set 
fi  in  r'^  represents  the  surrounding.  There  are  many  possible  boundary 
conditions  on  the  interface  between  the  slab  and  its  surroundings. 

1.  One  possibility  is  to  assume  that  the  slab  is  insulated  from  its 
surroundings.  This  condition  would  be  fulfilled  if  the  normal  derivative 
of  the  temperature  (outward  relative  to  the  set  ?.)  is  set  to  zero  on 

the  boundary  T  of 

2.  Another  possibility  is  to  suppose  that  the  surrounding  is  an  infinite 
sink,  i.e.,  its  temperature  is  not  affected  by  the  heat  diffusion  in 
the  slab.  This  state  would  be  Indicated  by  setting  the  temperature  of 
the  slab  on  the  boundary  T  to  a  constant,  which  may  be  assumed  to  be 


zero  by  translating  it. 


Many  other  possibilities  may  occur  such  as  a  combination  of  1  and 
2.  In  the  sequel,  condition  2  is  assumed  to  fix  the  ideas. 

The  smallness  of  e  means  physically  that  the  relative  diffusivity 
of  medium  1  is  small  with  respect  to  the  diffusivity  of  medium  0.  Letting 
£->■0,  e.g.,  in  Example  2.1,  signifies  that  medium  1  is  less  and  less 
conductive.  In  the  limit,  it  becomes  an  insulator.  This  situation  is 
symbolized  mathematically  by  the  normal  derivative  (outward  relative  to  2^) 

of  the  temperature  of  medium  0  going  to  zero  on  the  interface  S  of  the  two 

k.1 

media.  Consequently,  some  of  the  eigenvectors  of  ,  i.e., 

2 

do  reflect  this  behavior  as  indicated  by  their  weak  limits  (in  L  (:'))  given 
by  (2.32). 

From  the  viewpoint  of  medium  1,  medium  0  is  so  conductive  that  it 
may  be  considered  an  extension  of  the  surrounding  for  small  values  of  e.  If 
medium  0  is  insulated  from  (respectively  connected  to)  the  surrounding,  it 
becomes  an  insulator  (respectively  a  sink)  in  the  limit.  These  situations 
are  clearly  depicted  by  the  limits  of  some  of  the  eigenvectors,  i.e.,  {'f  } 

G  X 

in  Example  2.3  (respectively  Example  2.1). 

Now  consider  the  eigenvalue  problem  Example  2.6.  In  this  case,  the 
conductivities  of  both  media  are  of  the  same  order  of  magnitude.  However, 
the  convection  coefficients  are  of  0(1)  (respectively  0(e))  in  medium  0 
(respectively  medium  1),  i.e.,  although  the  heat  diffuses  in  the  slab  with 
comparable  rates,  the  internal  heat  exchange  with  the  surrounding  in  medium  0 
is  much  greater  than  in  medium  1  and  this  causes  stiffness.  Consequently, 
medium  0  is  a  better  heat  dissipator  than  medium  1. 

From  this  discussion,  it  seems  logical  to  expect  that  the 


eigenvalue-eigenvector  pairs  of  would  reflect  this  behavior  as  c-0. 


It  was  shown  that  the  eigenvalues  of  can  be  decomposed  into  two 
00  Ic 

groups  on  how  they  converge  as  e-^0. 

Exactly  as  previously  indicated,  medium  0  becomes  an  extension  of  the 
surrounding  as  e 0  and  eventually  a  sink  in  the  limit.  By  contrast, 
medium  1  loses  more  and  more  of  its  ability  to  dissipate  energy  as  e->-0, 
which  becomes  negligible  for  small  values  of  e.  This  is  clearly  demon¬ 
strated  by  the  limits  of  the  eigenvectors  given  by  (2.65a-b),  which  are 
completely  decoupled. 

Identical  interpretations  can  be  advanced  in  the  field  of 
electromagnetics,  provided  diffusivity,  sink,  insulator,  etc.,  are  replaced 
by  appropriate  terminology. 


5.3.  Numerical  Analysis  of  Parabolic  and  Hyperbolic  Boundary  Value  Problems 

In  this  section,  the  boundary  value  problems  of  Examples  3.6-3. 7  are 

are  revisited.  The  exact  solution  of  each  problem  is  not  available  for  the 

reasons  previously  discussed.  In  the  sequel,  the  zeroth  order  approximations 

obtained  in  the  aforementioned  examples  are  compared  with  the  finite- 

dimensional  approximations  of  these  boundary  value  problems. 

The  set  n=  (-1,1)  is  divided  into  N  equal  intervals  of  length 

h=^.  The  roof  functions  are  selected  as  a  basis  for  the  finite- 

N  h  i=l 

dimensional  approximation  of  Hp(n)  [2,20,39]. 

5.3.1.  Parabolic  problem 

It  is  straightforward  to  show  that  the  solution  of  the  boundary 


value  problem  (3.73)  can  be  approximated  by 


where  =  [c^  ...  (t  =  transpose)  satisfies 

M^c  +  k\  =  (5. 

e  e  e 

c^(0)  =  a^  (5. 

is  given  by  (2.89) 
is  given  by  (2.90) 


^i  "  i=  1.2....,N-1  (5. 

=  (a,^^),  i=  1,2,...,N-1.  (5. 

Remark  5.1;  Note  the  notation  change  from  Chapter  3,  i.e.,  =  g  on  H 

instead  of  7^(0)  = h  on  Q  because  h  designates  the  mesh  size.  ° 

Remark  5.2;  The  solution  of  (5.2)-(5.3)  for  t£(0,10)  is  obtained  by 
using  the  integration  routine  DGEAR  from  the  IMSL  library.  a 

Remark  5.3:  In  the  forthcoming  plots,  the  broken  lines  represent  y^(x,t) 
as  computed  in  Example  3.6  and  the  solid  lines  depict  y^(x,t).  □ 


For  all  computer  runs,  N  was  chosen  to  be  60.  Table  5.1 
summarizes  the  computer  runs.  It  is  noteworthy  to  mention  that  these  plots 
are  both  f inite~dimensional  approximations  of  y^  because  y^(x,t)  is  also 
approximated  by  a  finite  summation  (large  enough  to  obtain  a  smooth  plot!). 

It  is  evident  that  these  approximations  are  close.  For  small  e, 
they  almost  coincide  with  each  other.  Note  y^(x,t)  and  y^(x,t)  for  t  =  4, 
e*0.1  in  Figure  5.2a  are  at  steady  state.  The  temperature  distribution  on 
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TABLE  5.1.  PARABOLIC  PROBLEM 


Figure 

Plots  of  y 

e 

(x,t)  and  y^ 

(x,t)  for 

5.1a 

mam 

xe  (-1.1), 

e  =  0 . 1 

5.1b 

■B 

xe:  (-1,1) , 

e  =  0.001 

5.2a 

t  =  4, 

xe  (-1,1) , 

e  =  0.1 

5.2b 

t  =  4. 

x£  (-1.1). 

e  =  0.001 

5.3a 

t€  (0,10), 

X  =  -0. 5 , 

e  =  0.1 

5.3b 

te  (0.10), 

X  =  -0.5, 

e  =  0.001 

5.4a 

te  (0,10), 

X  =  0.5, 

e  =  0.1 

5.4b 

t€  (0,10), 

X  =  0. 1 , 

e  =  0.001 

is  a  straight  line,  i.e.,  because  it  is  due  to  a  point  source  on  the 
interface  x  *  0. 

The  temperature  at  x  =  0.5  as  a  function  of  time  for  e  =  0.001  is 

very  small.  Hence,  the  temperature  at  x  =  0.1  is  plotted  in  Figure  5.4b. 

Remark  5.4;  As  a  general  rule,  approximates  y^  pointwise  much  better  in 

the  interior  of  Q,  away  from  the  interface.  To  substantiate  this  claim, 

plots  5.3a-5,4b  are  provided.  In  plots  5.4b,  even  though  t  decreased  by  a 

0  h 

factor  of  100,  the  error  between  y^  and  y^  at  x  =  0.1  is  comparable  to  the 
error  at  x  =  0.5  for  e  =  0.1. 

5.3.2.  Hyperbolic  problem 

As  with  the  parabolic  problem,  the  solution  of  the  hyperbolic 


boundary  value  problem  (3.132)  can  be  approximated  by 


(5.6) 


Remark  5.5;  Cf.,  Remarks  5. 1-5.3.  □ 

For  all  computer  runs,  N  was  selected  to  be  60.  Table  5.2  sum¬ 
marizes  the  computer  runs.  Due  to  the  asymptotic  error  estimate  (3.11^)>one 

0  li 

would  expect  that  and  y^  would  not  be  as  "close"  as  in  parabolic  systems. 
Nevertheless,  the  two  approximations  "approach"  each  other  as  e-0.  These 
facts  are  clearly  illustrated  by  the  plots  in  Figures  5.5a-5.10b. 

Remark  5.6:  At  t=  1,  y^  is  growing  up  until  it  reaches  its  maximum  at  t  =  2. 
At  t  =  4,  it  attains  the  minimum  and  this  process  is  repeated  periodically 
every  four  units  of  time.  However,  behaves  differently  because  it  is 


Figure 


Plots 

of  y^(x,t)  and  y^(x,t)  for 

t  =  1. 

xe(-l,l).  E  =  0.1 

t  =  1. 

xe(-l,l),  £  =  0.001 

5.6a 

m 

,  xG  (-1 , 1)  ,  £  =  0.1 

5.6b 

t 

=  2 

,  xG(-l,l),  £  =  0.001 

5.7a 


5.7b 


5.8a 


5.8b 


5.9a 


5.9b 


5.10a 


5.10b 


;  =  4,  xG(-l,l),  £  =  0.1 


xG(-l,l).  £  =  0.0 


t  =  8,  xG(-l,l),  e  =  0.1 


t  =  8,  xG(-l,l).  £  =  0.001 


t€(0,10),  x  =  -0.5,  £  =  0.1 


t€(0,10),  x  =  -0.5,  £  =  0.001 


t€(0,10),  x  =  0.5,  £  =  0.1 


te(0,10),  x  =  0.1  £  =  0.001 


smitted  wave"  from  region  0  to  region  1.  Since  the  "velocity  of 


5. 


A  Parabolic  Control  Problem 


IbO 


In  this  section,  the  control  problem  of  Example  4.1  is  solved 
numerically.  Basically,  the  zeroth  order  approximation  y^  of  the  solution 
of  (4 . 58) -(4 . 61)  for  v  =  0  is  compared  with  the  controller  zeroth  order 

qC 

approximation  y^  of  C4.58)-(4.61)  for  v  given  (4.62).  •  The  various  results 
are  summarized  in  Table  5.3. 


TABLE  5.3.  PAPJ^OLIC  CONTROL  PROBLEM 


Plot 

Plot  of  y® 
e 

0*^ 

(x,t)  and  y^  (x,t)  for 

5.11a 

t  =  0.1, 

xe  (-1 , 1)  ,  e  =  0. 1 , 

p  =  0.1 

5.11b 

t  =  0.1, 

xG(-l,l),  £  =  0.001, 

p  =  0.1 

5.12a 

t  =  0.5, 

x€  (-1 , 1)  ,  £  =  0.1, 

p  =  0. 1 

5.12b 

t  =  0.5, 

xe(-l,l),  £  =  0.001, 

p  =  0.1 

5.13a 

t€  (0,1), 

X  =  -0 . 5 , 

p  =  0.1 

5.13b 

t€  (0,1), 

X  =  -0 . 5 , 

0  =  0.01 

5.14a 

te  (0,4), 

X  =  0.5,  £  =  0.1, 

p  =  0.1 

5.14b 

te  (0,4), 

x  =  0.5,  £  =  0.001, 

p  =  0.1 

Now  some  general  observations  are  in  order.  First,  if  no  control 


is  applied,  the  time  constants  associated  with  region  1  become  larger  and 
larger  as  z  decreases.  Hence  it  takes  longer  and  longer  for  the  state  to 
decay  to  zero.  This  fact  is  clearly  demonstrated  by  Figures  5.11b-5.13b. 


Second,  no  matter  how  small  e  is,  the  effectiveness  of  control  is 


not  appreciably  diminished.  This  is  readily  visible  in  Figures  5.12b, 
5.14b  and  may  be  substantiated  by  inspecting  the  Riccati  gains  given  by 
(4. 67), (4. 69) . 

Third,  the  control  is  regular  in  the  interior  of  Q,  away  from 
the  interface  as  shown  in  Figures  5.13a-5.14b. 

Finally,  the  effectiveness  of  the  control  is  enhanced  as  p 
decreases,  i.e.,  if  the  control  becomes  "cheaper,"  its  action  is  more 
effective  as  seen  in  Figures  5.13a-b.  This  is  not  a  pecularity  of  this 
example,  but  a  general  principle  of  control  theory. 

Remark  5.8;  The  effect  of  setting  y^j^  to  zero  on  R  produces  the  dip  in 
Figures  5.11a,  5.12a.  As  previously  indicated,  such  action  simplifies  the 
feedback  control  synthesis  and  induces  an  error  no  larger  than  0(v''T). 


5.5.  Concluding  Remarks 

In  this  chapter,  the  numerical  analysis  of  Examples  (3.6)-(3.7), 
(4.1)  is  undertaken.  The  conclusion  is  that  the  numerical  results  agree 
quite  well  with  the  theoretical  ones  obtained  in  Chapters  3  and  4. 

The  simplicity  of  the  examples  investigated  vV’ncea  1  s  .isf'ci  ts. 

For  example,  if  the  geometry  of  /,  is  more  complicated,  e.g.,  n-d imens i onal 
(n >  1)  and  polygonal,  the  direct  computation  of  even  the  limits  of  the  eigen 
vectors  is  very  elaborate.  Hence,  how  can  the  apnn'a.  h  t..  •  n  be 

extended?  It  was  indicated  in  Remark  5.7  that  a  direct  ipprsi,  h  osi  ig  i 


finite  element  method  is  hopeless.  Moreover,  verv  complex  intecr.jtion 
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routines  have  to  be  used  because  of  the  inherent  stiffness  of  the  problem 
at  hand. 


The  following  approach  seems  to  be  the  logical  alternative,  which 
consists  of  the  combination  of  the  direct  approach  and  the  approach  pursued 
in  this  thesis,  i.e.. 


^tep  1;  Use  a  finite  element  method  to  find  the  limits  of  the  eigenvectors 

Ic 

of  A  as  e->0.  In  so  doing,  the  limits  of  {y  },  ,  and  the  rate  of  conver- 

e  ®  e  k.=  l 


gence  of  { X 


£'k=l 


are  obtained  as  a  by-product  of  this  computation. 


Step  2;  Use  the  zeroth  order  approximations  derived  in  Chapter  3,  where 
Ic  Ic  Ic  k 

^k=l  ^k-1  replaced  by  those  computed  in  Step  1. 


In  the  first  step,  the  stiffness  of  the  problem  is  alleviated. 


In  the  second  step,  the  use  of  expensive  integration  routines  is  eliminated. 
However,  some  integration  is  still  required  to  get  the  desired  approximation. 


CHAPTER  6 


CONCLUSION 


6.1.  Concluding  Remarks 

In  this  report,  the  spectral  decomposition  of  some  stiff  partial 
differential  operators  is  undertaken.  One  class  of  such  operators  has 

coefficients  that  are  0(1), 0(e), ,0(e^)  in  . (r;^n  =  $  if  iii^j), 

whose  union  constitutes  the  open  connected  set  CCr’^.  It  is  found  that  the 
eigenvalues  of  these  operators  can  be  separated  into  p+1  groups,  depending 
on  how  they  converge  as  e-»-0.  Each  group  is  associated  with  a  region 
0^,  i=0,l,...,p.  The  eigenvectors  can  also  be  classified  accordingly. 

Their  convergence  as  c  -*■  0  is  intimately  related  to  the  order  of  the 
operators  and  the  location  of  the  coefficients  in  question.  A  general 
rule  has  emerged  out  of  the  present  investigation.  Suppose, without  loss  of 
generality, that  the  coefficients  in  question  occur  in  the  kth  derivatives 
(|k|_<  operator  order).  Then  the  higher  the  |k|,  the  weaker  is  the  convergence 
of  these  eigenvectors  as  e  ^  0.  This  conclusion  is  not  surprising  because 
in  this  case  the  eigenvectors  are  more  regular, in  general  [30] .  The  details 
for  some  typical  cases  are  summarized  in  Table  2.1. 

In  this  thesis,  the  value  of  the  parameter  e  is  assumed  to  be 
small.  Consequently,  the  behavior  of  the  spectrum  of  stiff  operators  is 
analyzed  as  e  0.  However,  similar  results  can  be  stated  as  e  ->•  +  “,  as 
evidenced  by  Corollary  2.1. 


In  this  thesis,  it  is  also  shown  that  approximations  of  solutions 
of  stiff  boundary  value  problems  (including  some  control  problems)  can  be 
derived  using  the  weak  limits  of  the  eigenvectors  of  stiff  operators  as 
£  ->■  0.  These  approximations  are  "readily"  computable.  However,  in  general, 
they  are  not  as  regular  as  the  exact  solutions,  which  may  not  be  accessible 
at  all. 

6.2.  Topics  for  Future  Research 

Many  unsettled  issues  directly  related  to  this  thesis  merit  farther 
research,  some  of  which  are: 

1.  The  analyticity  of  (Cf.  Remark  2.26). 

£  iC“  X 

2.  When  is  not  connected,  e.g.,  Slg  =  Qgg  U  12^^  (Cf.  Remark  2.10), 

Ic  Ic  Ic  Ic 

the  limits  of  »  i.e.,  ^k=l  decomposed  further  into 

k  k  ®  k  Ic  ® 

{Uoo»'»'o0^k“l  ^^01’'^01^k-l’  pair  associated  with  a  subset  of  .Ig,  as 

indicated  by  the  subscripts.  It  is  not  clear  if  there  exists  one  or  two 
eigenvalue-eigenvector  pairs  for  c  >  0,  which  correspond  to  each  of  the 
above  pairs. 

Finally,  it  would  be  worthwhile  tO  examine  if  the  present 
methodology  can  be  used  to  approximate  the  solutions  of  the  following 
problems : 

3.  Semilinear  boundary  value  problems  [4], 

4.  Unilateral  problems  [4,23], 

5.  Inverse  problems  [6], 

6 .  Games  [3,35]  , 


when  they  involve  stiff  operators. 


APPENDIX 


In  this  appendix,  some  of  the  mathematical  tools  needed  in  this 
thesis  are  discussed  and  the  inquiring  reader  is  referred  to  the  appro¬ 
priate  references. 

The  following  subjects  are  very  briefly  reviewed. 

1.  Definition  of  0(0  and  o(0 

2.  Weak  convergence 

3.  Distributional  derivatives  and  functional  spaces 

4.  Definition  of  the  space  £ (X;Y) 

1.  Def •’’nition  of  0(0  ar.d  o(0 

Let  fg(t)  and  fj^(£)  be  real,  positive,  continuous  functions  of 

in  0  <  £  <  such  that  lim  f.(e),  i  =  0,l  exist. 

■  ^  e^O  ^ 

1.1.  Definition  of  0(0 

fg(£)  =  0(fj^(£))  if  there  exists  a  constant  C  such  that 
fQ(£;)  <  Cfj^(£)  for  0<£  <eQ. 

1.2.  Definition  of  o(0 

f„{£) 

£„(,;)  -  0(£^(c))  if  -  0 

For  more  details,  consult  [11]. 

2.  Weak  Convergence 

Let  H  be  a  real  Hilbert  space  with  scalar  product  (x,y)^. 


sequence  -x  •  converges  weakly  to  an  element  x£H  if 
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(Xn.y)^ (x.y)H' 

For  further  inquiry,  consult  [42]. 


3.  Distributional  Derivatives  and  Function  Spaces 

Let  X  =  {Xj^  jX^  j  •  • .  ,x^}  denote  the  space  variable;  x  ranges  over 

an  open  set  with  boundary  F;  t  denotes  time,  t£(0,T),  T  <  “.  Let 

C  (Q)  =  space  of  k-tiraes  continuously  differentiable  functions  on  7.,  kCiN 
Ic 

C^(Q)  =  space  of  k-times  differentiable  functions  in  7,  with  compact  support 


in  Q . 

P(7)  =  Cq(7). 

V*  (Ti)  =  dual  space  of  0(7),  i.e.,  the  space  of  distributions  on  7. 

L^(r.)  =  space  of  functions,  which  are  square  integrable  on  7,  which  is  a 

2 

Hilbert  space  (L  (7)  is  identified  with  its  dual) . 


3.1.  Distributional  derivatives 


Let 


p  =  {p^,p,,...,p^},  !p!  =  p^ 


+  ?„  +  •••+? 


p  Pi  P2  Pn 

=  D,  D„  .  .  .  D  ,  D.  =  -r^ 
12  n  1  jx . 


f€L"(7) 


is  said  to  be  the  qth-distributional  derivative  of  f  if 


(D'’f,r')  =  (-1) <f ,  Vv-e0(7) 


where  designates  duality  between  0(7)  and  O' (7) . 


One  may  define  the  Sobolev  space  of  order  m  as 


r  ^  *  2  ^  »  2  , 


rxfl)  =  :v'e  L"(Q),  d\£L  (fi),  Vq,  ]  q  |  <  ml 


which  is  a  Hilbert  space  if  endowed  with  the  scalar  product 


.m 


H  (2)  I q 1  <  m 


For  any  v^SH^CZ),  one  can  define  uniquely  its  traces  on  the  boundary,  i.e.. 


(  -^1 

V  «  Tt"  - - - 

Ir 


ra-1 

3 _ ^ 


r 


3v 


m-1 


k  I 

■3  s? 


where  s?  ' - »  { — j-  »  0  :£  k  <  m-1}  is  a  continuous  linear  surjective  map 

3v^  I  r 


of  onto  •  Using  the  trace  map,*one  can  define  several 

k=0 


subspaces  of  H™(r:)  such  as 


H"'(a;r  )  =  :  v'e  h"'(Z),  ~ 


=  0 ,  FqC  r ,  0  <  k  <  m- 1 } 


‘0 


hJcz)  =  H“(r.;r) 


If  time  is  involved,  many  Hilbert  spaces  can  be  defined  in  a  similar  way, 
2, 


e.g.,  L  (0,T;V)  =  space  of  functions  defined  on  (0,T)  with  values  in  a 
space  V  such  that 


0 


where  V  may  be  any  of  the  above  Sobolev  spaces  or  their  duals. 
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If  V  =  (ro  ,  then 

L“(0,T;V)  =  L^(O.T;L^(C))  =  (0,T)) 

A  systematic  study  of  these  spaces, as  well  as  many  of  their  subspaces,  is 
found  in  [2,30]. 


4.  Definition  of  £  (X;Y) 

jC(X;Y)  is  the  vector  space  of  continuous  linear  operators  from  X 
to  Y.  The  norm  of  an  operator  AGT(X;Y)  is  defined  by 


llAll 


sup 

xex 


IIAxIIy 

l|x||^ 


If  Y  is  a  Hilbert  space  and  X  is  a  pre-Hilbert  space,  the  space  (X;Y) 
is  a  Banach  space.  See  [2,15,42]  for  more  details. 
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